Post-Newtonian, Quasi-Circular Binary Inspirals in Quadratic Modified Gravity 
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We consider a general class of quantum gravity-inspired, modified gravity theories, where the 
Einstein-Hilbert action is extended through the addition of all terms quadratic in the curvature 
tensor coupled to scalar fields with standard kinetic energy. This class of theories includes Einstein- 
Dilaton-Gauss-Bonnet and Chern-Simons modified gravity as special cases. We analytically derive 
and solve the coupled field equations in the post-Newtonian approximation, assuming a comparable- 
mass, spinning black hole binary source in a quasi-circular, weak-field/slow-motion orbit. We find 
that a naive subtraction of divergent piece associated with the point-particle approximation is ill- 
suited to represent compact objects in these theories. Instead, we model them by appropriate 
eiTective sources built so that known strong-field solutions are reproduced in the far-field limit. In 
doing so, we prove that black holes in Einstein-Dilaton-Gauss- Bonnet and Chern-Simons theory can 
have hair, while neutron stars have no scalar monopole charge, in diametrical opposition to results 
in scalar-tensor theories. We then employ techniques similar to the direct integration of the relaxed 
Einstein equations to obtain analytic expressions for the scalar field, metric perturbation, and the 
associated gravitational wave luminosity measured at infinity. We find that scalar field emission 
mainly dominates the energy fiux budget, sourcing electric- type (even-parity) dipole scalar radiation 
and magnetic-type (odd-parity) quadrupole scalar radiation, correcting the General Relativistic 
prediction at relative — IPN and 2PN orders. Such modifications lead to corrections in the emitted 
gravitational waves that can be mapped to the parameterized post-Einsteinian framework. Such 
modifications could be strongly constrained with gravitational wave observations. 

PACS numbers: 04.30.-w,04.50.Kd,04.25.-g,04.25.Nx 



I. INTRODUCTION 

The validity of Einstein's theory in the strong-gravity 
regime will soon be put to the most stringent tests yet, 
through the observation of gravitational waves (GWs) 
from compact object binary inspirals [J-Q- Such waves 
carry detailed information about their source and the 
underlying gravitational theory in play. This informa- 
tion is primarily encoded in the evolution of the GW 
frequency, which in turn depends directly on the rate of 
energy transport away from the binary In general 
relativity (GR), this transport is performed exclusively 
by GWs. In modified gravity theories, however, addi- 
tional (scalar, vectorial or tensorial) degrees of freedom 
can also carry energy and angular momentum away as 
they propagate. 

Calculating how gravitational waves are corrected in 
modified gravity theories can be a gargantuan task as 
the modification can increase the number of propagating 
degrees of freedom and the non-linearity of the equations 
that control their propagation. For example, the amount 
of energy-momentum transported away from a binary 
system must be computed both from the GWs excited 
by the corresponding sources, as well as any additional 
waves associated with extra degrees of freedom Q. The 
sources that drive such waves can depend both on deriva- 
tives of the metric perturbation and the extra degrees of 
freedom, which, in turn, arc specified by the solution to 
their own equations of motion. The situation worsens if 



these are non-linearly coupled, e.g. a scalar field equation 
of motion that depends on the metric tensor, whose evo- 
lution in turn depends on derivatives of the scalar field. 

Such calculations, however, are feasible if one treats 
any GR deviations as small deformations , which can 
be formalized through the small- coupling approximation, 
a common technique in perturbation theory to isolate 
physically relevant solutions in higher-derivative theo- 
ries • This is a reasonable approximation given that 
GR has passed a large number of tests, albeit in the 
weak-gravity regime. Even in the GW regime, signals 
will slowly transition from sampling weak fields to mod- 
erately strong fields during a full binary inspiral. The 
strongest GW events will not be able to sample anywhere 
close to the Plank regime, where one would expect com- 
pletely new physics. The largest gravitational fields ex- 
perienced by binaries occur when these merge, and even 
then, the metric curvature cannot exceed m~^, where m 
is the total mass of the binary. Earth-based detectors, 
such as LIGO [l^, VIRGO [ij and LCGT Jp, and fu- 
ture space-borne detectors, such as LISA [ij], will only 
be able to sample gravitational fields up to this strength. 

Of the plethora of modified gravity theories, we choose 
to focus on a general class that is characterized by the 
addition of quadratic curvature invariants to the ac- 
tion, coupled to scalar fields with standard kinetic terms 
(see e.g. Eq. ([T|). Such theories are motivated from 
loop quantum gravity [HI and heterotic string the- 
ory [16[, arising generically upon four-dimensional com- 
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pactification in the low-energy limit. Disjoint sub-classes 
of quadratic theories reduce to Einstein-Dilaton-Gauss- 
Bonnet (EDGE) theory [13, [lland Dynamical Chern- 
Simons (CS) modified gravit70, ■ 

From a phenomcnological standpoint, such quadratic 
gravity theories are also interesting as straw-men to 
study small deviations from GR. This is because the 
new quadratic terms are always small relative to the 
Einstein-Hilbert term when considering merging binaries. 
In such systems, the minimum radius of curvature is al- 
ways larger than the new scale introduced by the scalar 
fields. If this were not the case, astrophysical observa- 
tions would already have constrained quadratic gravity 
deviations. 

Quadratic gravity introduces an equation of motion 
for the scalar field and modifies the metric field equa- 
tions. The former is a driven wave equation, whose 
sources are quadratic curvature invariants. The latter 
contains new terms that depend on the product of the 
scalar field and its derivatives with the Riemann tensor, 
Ricci tensor, Ricci scalar and their derivatives. As such, 
one might worry that higher derivative terms in the field 
equations could render the system unstable. One must 
remember, however, that the action is a truncation (at 
quadratic order in the present case) of an effective the- 
ory derived by integrating out heavy degrees of freedom 
contained in a more complete theory. Since we truncate 
the effective action, its validity is limited only to leading- 
order in the coupling parameters. Accounting for higher- 
order terms in the coupling would require the inclusion 
of higher-order terms (cubic, quartic, etc.) in the ac- 
tion [8|. Therefore, the modified field equations should 
not be considered as an exact system, but rather as an 
effective one. 

Given the above and using the small- eoupling approx- 
imation, the field equations become driven differential 
equations for the metric deformation and the scalar field. 
The source of the latter depends only on derivatives of 
the GR metric perturbation, while the source of the 
former depends both on the GR metric perturbation 
and the scalar field. We solve these equations in the 
post-Newtonian (PN) limit, where in particular we con- 
sider comparable-mass, spinning black hole (BH) bina- 
ries (electromagnetically uncharged), spiraling in a quasi- 
circular orbit. This forces the driven differential equa- 
tions into driven wave equations, which can be studied 
with PN techniques [2l| - |27| and then solved via retarded 
Green function methods. 

A complication arises when attempting to solve these 
equations, as one must choose a prescription to de- 
scribe BHs and neutron stars (NSs). In standard PN 
theory and up to a certain high PN order, one can 
choose a point-particle prescription, essentially because 
the exterior gravitational field of a compact object is 
the same as that induced by a point-particle. In mod- 
ified quadratic gravity, however, both non-spinning [28j 
and spinning [6|, strong- field BH solutions differ from 
that generated by simple point particles with a mass- 



monopole and a current-dipole moment; BHs in these 
theories have additional scalar multipole moments. (See 
Refs. m, 1131 for similar discussions on NSs in CS grav- 
ity.) One can take these effects into account by con- 
structing an effective point-particle source that repro- 
duces known, strong-field solutions to leading order in 
the weak-field region, sufficiently far away from the com- 
pact objects. With this effective point-particle prescrip- 
tion, we can then evaluate the source of the driven wave 
equations and analytically solve them to find the radia- 
tive part of the scalar field and metric perturbation. 



Executive Summary of Results 

Given the length of this paper, let us summarize the 
main results. We have devised a framework in the small- 
coupling approximation to solve for compact binary in- 
spirals in modified quadratic gravity theories. One of the 
key ingredients in this framework is the calculation of ef- 
fective source terms that allow us to use the point-particle 
approximation even for theories where such approxima- 
tion is not valid. We applied this to modified quadratic 
gravity to find that both NSs and BHs have scalar hair, 
which leads to dipolar emission. EDGB and CS grav- 
ity are exceptions, where although BHs retain scalar 
monopole and dipole charge, respectively, NSs shed the 
scalar monopole charge. Therefore, BHs in EDGB gener- 
ically contains dipolar GW emission, while CS gravity 
leads to modified quadrupolar emission. 

The presence of scalar monopole and dipole hair, and 
in particular the fiux of energy-momentum carried by 
this hair, leads to a modification in the rate of change 
of the binary's binding energy. The even-parity sector 
of the theory leads to scalar hair, which modifies the en- 
ergy flux at — IPN order relative to the GR quadrupole 
fiux. Of course, such a modification is proportional to 
the coupling parameter of the theory, which is assumed 
small. The odd-parity sector leads to dipole hair for spin- 
ning BH binaries, which modifies the energy flux at 2PN 
relative order. If the BH binary components are non- 
spinning, they have no dipole hair but the binary orbital 
interaction generates a modification in the energy flux 
that enters at relative 7PN order. Figure [T] shows the 
energy flux carried by the even-parity scalar field (long 
dashed line), odd-parity scalar field (dot-dashed for spin- 
ning binaries and short dashed line for non-spinning bi- 
naries), and the GR quadrupole flux (solid line) as a func- 
tion of orbital velocity. Observe that when one assumes 
that BHs are non-spinning, the scalar emission is greatly 
suppressed. 

These energy flux corrections translate into changes to 
the waveform observables. We explicitly calculate these 
and map them to the parametrized post-Einsteinian 
(ppE) framework [3l|, [32|. Using the results of Cor- 
nish et al. [32| we estimate that GW observations could 
constrain the new length scale introduced in quadratic 
gravity (related to the coupling constants of the the- 
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FIG. 1: Comparison of the energy flux carried by scalar fields 
of even-parity (dashed red), odd-parity and sourced by spin- 
ning BHs (blue dot-dashed) and odd-parity and sourced by 
non-spinning BHs (short dashed) relative to the GR pre- 
diction (solid black) as a function of orbital velocity. We 
here consider a quasi-circular, BH inspiral with (mi, 7712) = 
(8,2O)M0, normalized spins S\ = \Sl\/ml = -S2 = 
— \S2\/m2 perpendicular to the orbital plane, 15*^1 — m\ and 
coupling constants (^3 = 6.25 x 10~^ = ("4. 

ory) to roughly the BH horizon scale. With a typi- 
cal Ad. LIGO stellar-mass BH inspiral observation, one 
should be able to constrain the even-parity sector to 
roughly 0(10) km. With a typical LISA extreme-mass 
ratio inspiral (EMRI) observation, one should be able to 
constrain the odd-parity sector to roughly 0(100) km. 
Such projected constraints are much stronger than cur- 
rent Solar System bounds [13, l33l - [35j . 

This paper is organized as follows: Section HIl describes 
the action that will be considered in this paper and re- 
views the associated modified field equations and the 
scalar field equation of motion. Section IIIII expands the 
field equations in the small-deformation approximation. 
Sections ll VI and fVl study the scalar field and metric defor- 
mation evolution, analytically solving the modified field 
equations. Section |VI] computes the energy flux carried 
by the scalar field and the metric deformation. Sec- 
tion |VlT] considers the impact that such fluxes would have 
on gravitational waveform phase. Section IVlIII concludes 
and points to future research. 

We have deferred many details of the computational 
techniques to the appendices. Appendix [Xlshows the NSs 
in EDGB theory have no scalar monopole charge. Ap- 
pendix |B] discusses speciflc integration techniques. Ap- 
pendix [C] estimates the order of the metric correction 
from the regularized contribution for non-spinning BHs 
in the odd-parity sector of the modified theory. Ap- 
pendix |D] discusses particular integrals that appear when 
solving the field equations. 

Henceforth, we follow mostly the conventions of Mis- 
ner, Thorne and Wheeler Greek letters stand for 
spacetime indices; Latin letters in the middle of the al- 
phabet i, j, . . ., stand for spatial indices only. Parenthe- 
sis, square brackets and angled brackets in index lists de- 



note symmetrization, antisymmetrization and the sym- 
metric and trace free (STF) operator, respectively. Cap- 
ital Latin letters usually refer to a multi-index, such as 
rj.Q — x'^jk... ^ where x*^*^" = x^x^x'^ .... Partial deriva- 
tives are denoted with diA = A^i = dA/dx"^, while co- 
variant derivatives are denoted with the nabla ViA, for 
any quantity A. Deformations arc labeled with the order- 
counting parameter Finally, we use geometric units, 
where G = c = 1. except when denoting the order of 
certain terms in the PN approximation. Throughout, we 
performed analytic calculations with the xTensor pack- 
age for Mathematica [H, [13 . 

II. MODIFIED GRAVITY THEORIES 

In this Section, we introduce the class of modified grav- 
ity theories that we study, by writing down its action and 
equations of motion. We then proceed to define the small 
deformation approximation more precisely. 

A. ABC of quadratic gravity 

Consider the following 4-dimcnsional effective action: 
S = J d*x^/^{Ki^ + al/l(^?)i?2 + a2/2(^?)i?p.i?^" 

- f [v^i^v^T? + 2v{d)] + /:,„at| . (1) 

Here, g stands for the determinant of the metric .g^jy. 
Rnv, RfiuSa and * R^vSa are the Ricci scalar and tensor, 
the Riemann tensor and its dual [ssf . respectively, with 
the latter defined as^ *R^^5a = {l/2)e5a°''^ Rf'^afs and 
with et^"^'^ the Levi-Civita tensor. The quantity £mat is 
the external matter Lagrangian, is a field, {ai,f3) are 
coupling constants and k = (IBtt)"^. This action con- 
tains all possible quadratic, algebraic curvature scalars 
with running (i.e. non-constant) couplings, where we as- 
sumed that all quadratic terms are coupled to the same 
field. All other quadratic curvature terms are linearly 
dependent, such as the Weyl tensor squared. 

The theory defined by the action above is different 
from f{R) theories on several counts. First, f{R) the- 
ories depend only on the Ricci scalar, while the action 
above depends on the Ricci tensor, the Riemann tensor 
and a dynamical field 1?. Second, f{R) theories are usu- 
ally treated as exact, while the action presented above is 
an effective theory, truncated to quadratic order in the 
Riemann tensor. The consequence of this is insisting on 
the use of order-reduction in the field equations, where 



^ This definition is correct, in agreement with [38l |. and fixing an 
inconsequential typo in [3911 . 
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we treat all quantities that depend on at perturbatively. 
Such order reduction then leads to the absence of addi- 
tional polarization modes [H, li^l , such as the longitudi- 
nal scalar mode that arises in f{R) theories. 

The field equations of dynamical quadratic gravity can 
be obtained by varying the action with respect to all 
fields. For simplicity, we restrict attention to coupling 
functions fi{'&) that admit the Taylor expansion fi{^) = 
/^(0) + //(0)i? + 0{i}^) about smaU i}, where /i(0) and 
f^(0) are constants, and we assume that the asymptotic 
value of 9 at spatial infinity vanishes. Let us further 
reabsorb /i(0) into the coupling constants af''^ = aifi{0) 

and /j'(0) into the constants a[^^ = aif-{0). Equation ^ 
then becomes 5* = S'gr 4- 5*0 + 5*1 : 



'5'gr — 
Si = 



J d^X^{K.R + Crn.t} , (2) 

+ af^R^^saR'^"'^} , (3) 
J d^xy/^^^ai^^-dR^ + ai^'^dRf.^R'"' 



2 



(4) 



where clearly is the Einstein-Hilbert plus matter ac- 
tion. Notice that Sq defines a GR correction that is de- 
coupled from 0. The term proportional to a][''^ can not af- 
fect the classical field equations since it is topological, i.e. 
the second Chern form, so we have omitted it. Similarly, 



.,(0) 



variant, (al , a 



(0) ^,(0) „.(o) 



a. 



) = (1, —4, l)aGB, then these will 



not modify the field equations. On the other hand, S'l 
defines a modification to GR with a direct (non-minimal) 
scalar field coupling, such that as the field goes to zero, 
the modified theory reduces to GR. We here restrict at- 
tention to the case aj*^^ = 0. From this point forward, we 
will drop the superscript from a[^\ 

The action above defines a class of modified gravity 
theories that contains well-known GR extensions. For ex- 
ample, when a4 = —jUcs and all other ai = 0, quadratic 
gravity reduces to dynamical CS gravity, where ctcs is 
the CS coupling parameter (see e.g. [s^). Alternatively, 
when a4 = 0, while (01,02,0:3) = (1, -4, l)aEDGB, 
quadratic gravity reduces to Einstein-Dilaton-Gauss- 
Bonnet theory (see e.g. [3). Both of these theories are 
motivated from fundamental physics; they unavoidably 
arise as low-energy expansions of heterotic string the- 
ory [4l| - |43 |. Dynamical CS gravity also arises in loop 
quantum gravity when the Barbero-Immirzi parameter 
is promoted to a field in the presence of fermions j45l - l47| . 

Variation of the action with respect to the metric yields 



the modified field equations: 



+^7^(1) + ^i(i) + ^ 7(1) + ^r(i) 

r\i r\i r\i Kj 

where we have defined the short-hands^ 
n^°J =2RR^, - ^g^,R^ ~ 2V^,i? + 2g^,ni? , 

^^4°' ^D-R^iy + 2R^Siy(yR^" ^ -^9 fii^ R^"^ RSa 

+ ^g^,DR-\/^,R, 
J(0) =8R^^R,,sua - 2g^,R^^Rsa + ADR^, 



(5) 



t(i) 



-t(i) 



(6a) 

(6b) 

(6c) 
(6d) 

(6e) 

8 {V^d) (V(^i?,)5 - VsR^,,) + 4i?^iv<.V*"i? , 

(6f) 



— 2RR^^ + -g^,^R — 2V^^i?, 

-A{V(^d)V,)R-2RV^,^ 
+ [2RUd + 4(V^^9)V5i?] , 
-(V(^i?)V,)i?-2V*W(^i?,)i- 
+ 2V'^i9V5i?^, + i?„,ni? - 2i?i-(^V*V,)^? 
+ g^, {V^-dVsR + R^^Vsa^) , 



/CW = - 4 (V^i?) esa^i^V'^R.f + 4{Vsa^rRi^'.f ' 

(6g) 

where is the covariant derivative, V^ii/ = V^Vi/, and 
n = V^V is the d'Alembertian operator. The •& field's 
stress-energy tensor is 



(V^i9)(V,i9) - ig/.- (ViW^i? - 2V{'d)) 



(7) 

Variation of the action with respect to d yields the d 
equation of motion: 

dij 

— OL^RfivScyR^^^'^ — OLiR^iySa R^^^'^ ■ (8) 

Notice that when the spacetime is curved by some mass 
distribution, the right-hand side will be proportional to 
density squared. 

The parity of the field i? can be inferred from its equa- 
tion of motion. Since terms of the form R? are even- 
parity, while terms of the form R^^sa* R'^'^^'^ are odd- 
parity, the field is of mixed parity. Note however that 
the even and odd-parity couplings tend to have different 



^ This corrects an error in Eq. (5b) of [28 
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origins from an underlying theory. In this paper we will 
consider the even and odd-parity cases separately. 

The inclusion of dynamics for the field in the ac- 
tion guarantees that the field equations are covariantly 
conserved without having to include any additional con- 
straints, i.e. the covariant divergence of Eq. ^ identi- 
cally vanishes, upon imposition of Eq. This is a 
consequence of the action being diffeomorphism invari- 
ant. Such invariance is in contrast to the preferred-frame 
effects present in a non-dynamical theory [l9| , i.e. in the 
theory defined by the action in Eq. (j?]) but with (5 = 0. 
In the latter, the field must be prescribed a priori. 
Moreover, the theory requires the existence of an addi- 
tional constraint (the right-hand side of ([8]) to vanish), 
which is an unphysical consequence of treating ■& as prior 
structure [48ll49l|. 

Before proceeding, let us further discuss the scalar field 
potential V{'d). This potential allows us to introduce 
additional couplings, such as a mass term, to drive the 
evolution in Eq. (|5]). However, there are reasons one 
might restrict such a potential. If the mass is much 
larger than the inverse length scale of the system that 
we concern, the effect of such a field on the dynamics of 
binaries is strongly suppressed. To the contrary, if the 
mass is much smaller, the presence of mass does not give 
any significant effects. Therefore we cannot expect to 
observe the effects of a finite mass without fine tuning. 
No mass term may appear in a theory with a shift sym- 
metry, which is invariance under d ^ d + const. Such 
theories are common in 4D, low-energy, effective string 
theories HI El H HO, Uli , such as dynamical CS and 
EDGB. For these reasons, and because the assumption 
makes the resulting equations analytically tractable, we 
will henceforth assume V{^) = 0. 



B. Small deformations 

The "unreasonable" accuracy of GR to explain all ex- 
perimental data to date suggests that it is an excellent 
approximation to nature in situations where the gravita- 
tional field is very weak and velocities are very small rela- 
tive to the speed of light. GW detectors will be sensitive 
to events in situations where the field is stronger than 
ever previously sampled. This, however, does not im- 
ply that GWs will ever sample the Planck/string regime, 
where one could expect large deviations from GR. 

We will here be interested in binary compact object 
coalescences up until the binary reaches the innermost 
stable circular orbit (ISCO). Even during merger, the 
largest curvature that GWs will sample will be limited to 
the scale determined by the horizon sizes, proportional 
to m~^. Such scales are far removed from high-energy 
ones, like the electroweak one, as GW detectors will not 
be sensitive to mergers of compact objects with masses 
below a solar mass. Even then, however, GWs can and 
will probe the strong field, which has not been tested 



before. One is then justified in modeling GWs that may 
contain deviations from GR as small deformations. 

The small deformation scheme is also appealing for 
theoretical reasons. As mentioned earlier, the theories 
we consider are effective, valid only up to the truncation 
order. There are higher-order terms that we have here 
neglected in the action, such as cubic and quartic cur- 
vature combinations. Thus, one should not treat these 
theories as exact nor insist on solving the equations of 
motion to higher orders in a^. If this is desired, then 
higher-order curvature terms should also be included in 
the action. 

One might be worried that such effective theories are 
unstable, since they lead to field equations with deriva- 
tives higher than second order. Such derivatives could 
lead to instabilities or ghost modes if the Hamiltonian is 
not bounded from below. Linearization in the coupling 
parameter, however, has the effect of recasting the field 
equations in Einstein form with an effective stress-energy 
tensor that depends on the GR solution, thus stabiliz- 
ing the differential equations Q. Linearization removes 
modes besides the two that arise in GR [39, 40]. 

Small deformations can be treated similarly to how one 
models BH perturbations. That is, we expand the metric 
as 

9^.^9^^ + <^^f.. + 0{^), (9) 

where the GR superscript is to remind us that this quan- 
tity is a GR solution, while is a metric deformation 
away from GR. The order-counting parameter <; is kept 
around only for book-keeping purposes and is to be set 
to unity in the end. 

Applying such an expansion to Eq. ([5]), one finds 

(3Dd = -a,S{RlJ + 0{<;), (10) 

where 5(i?gf^) stands for all source terms evaluated on 
the GR background g^^. The solution to this equation 
will obviously scale as oc ai//3. Applying the decom- 
position and expansion of Eq. ^ to Eq. ^ in vacuum, 
one finds 

G..[f)..] = --C,A^,9^^] + ^At¥] , (11) 

where the 0{<p) terms automatically vanish, as sat- 
isfies the Einstein equations, and we have grouped mod- 
ifications into the tensor C^^. This tensor and T^i^ are 
to be evaluated on the GR metric and act as sources 
for the metric deformation. Notice that, as a differential 
operator acting on f)^^, the principal part of these dif- 
ferential equations continues to be strongly hyperbolic, 
as it is still given by the G^^ differential operator, with 

the higher derivatives in C^^ and the T/ijy acting as 
sources. Given this, the metric deformation is propor- 
tional to = af/{l3K), which is our actual perturbation 
parameter. 
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Proper perturbation or deformation parameters should 
be dimensionless, but the are dimensional. The dimen- 
sions of a and f3, of course, depend on the choice of di- 
mensions for the scalar field. We here take the viewpoint 
that 1? is dimensionless, which then forces (3 to be dimen- 
sionless as well as k, and a to have dimensions of length 
squared. Then, the deformation parameter ^ has units of 
length to the fourth power, which is why we define the 
dimensionless 



(12) 



as our proper deformation parameter. One could choose 
different units for the scalar field, but in all cases one 
arrives at the conclusion that Q is the proper deformation 
parameter 



III. EXPANSION OF THE FIELD EQUATIONS 

Let us decompose the GR metric tensor into a flat 
background plus a metric perturbation: 



h„ 



(13) 



We emphasize here that throughout this paper, h^^, de- 
notes the metric perturbation in GR while is the 
metric deformation away from GR. 

In expanding the modified field equations, we will also 
find it useful to define the standard trace-reversed metric 
perturbation in GR as 



(14) 



In particular, notice that when the background is flat 
htiy = h^i, - ^hr]f,i, and /i^^ = hf^i, - ^hrifj,„ to hnear 
order in GR. We also define the deformed trace reversed 



metric perturbation as 



(15) 



The harmonic gauge condition reduces to /i^'^.j/ = and 
i)'^'^ ,1, = 0. Throughout this paper, we only study the GR 
deformation up to ©(a^//?) for ■& and 0(0) for fj^j,. 



A. Scalar field 

The evolution equation for the scalar field at leading 
order in the metric perturbation becomes 

2a 



(16) 



with relative remainders of 0{h). Here, e^'^^'^ is the 
Levi-Civita symbol with convention e^-*^^^ = +1 in an or- 
thonormal, positively oriented frame, and we have used 
the harmonic gauge condition. 



B. Metric perturbation 

Let us now perturb the metric field equations [Eq. ([5|)] 
about ^ 0. The deformed metric wave equation at 
linear order in becomes 



1 C^mat _ It^Ii?) 

2 /il^ rj /II/ ' 



(17) 



where the tensors on the right-hand side are given by 

- -4(/ip["'''l.p.-77,,.n^/ipI"^''l.) , (18) 

+7]^,D^hp^'''PK , (19) 



= ^{-a„h^^^,pf -a^hp^p,^], - h^'^'^K,.) , (20) 

?iW = -8/ip["'''l,(p^^,)-4/ip[^'''l,,?^p, 

+4?7/.. {2hp^^-PKs^'^ + hp^^'P^a^^^ , (21) 
fW ^ -2/ip[-'''l,(,,^,,)-2(/i*[,,(,]^)''-/i''[,,(,]^)*)7?,5 

-2 (V[m):P]<5'' + hP[p,(^^]^^s) 

-2{h\p,^,f^,^^s-hP[p.^p]'d,,^s) 

+r/p, [2h;"'PKs^-' + (/i^f^'^ip - V'"-'!") ^.as} 



+D^-d ( ,i,)s — -Hrfh^i, — -h,py j , 



(22) 



(23) 



-2^'^e5<,^p/i'^[„'"x,] + (Ai o y) , (24) 

where is the d'Alembertian of flat spacetime, h = 
/ip^, and Tpjy is given as 



(25) 



The quantity (5r^"J" stands for the perturbation to the 
energy-momentum tensor for matter. Even when dealing 
with BHs, STj^^* ^ because we treat BHs as distribu- 
tional point particles and their trajectories are generi- 
cally modifled at 0{<^). However, in this paper we con- 
centrate on the dissipative sector of the theory only, and 
not on modifications to the shape of the orbits (conser- 
vative dynamics) . The latter does modify the GW phase 
evolution [1,[2^, as we discuss in Sec. I Villi 
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The evolution equation for the metric perturbation 
takes on the same form (a sourced wave equation) as 
that for the scalar field. The source terms in both of 
these equations depend on the GR metric perturbation, 
which we here assume to be that of a compact binary 
quasi-circular inspiral in the PN approximation, i.e. mov- 
ing at small velocities relative to the speed of light and 
producing weak gravitational fields. We provide explicit 
expressions for the GR metric perturbation in the subse- 
quent subsection. 



C. Post-Newtonian metric and trajectories 

In this subsection, we provide explicit expressions for 
the linear metric perturbation in GR that we use to eval- 
uate all source terms. We are here interested in a binary 
system, composed of two compact objects with masses 
nil and m2 and initially separated by a distance ri2 = b. 
The objects' trajectories can be parameterized via 



X2 = X2 



TO2 



b [cosojt, sinwt, 0] , 



(26) 



6 [cos wt, sincji, 0] , (27) 

m 



where m = mi + m2 is the total mass and where we have 
assumed they are located on the x-y plane. Throughout 
this paper, vectors are sometimes denoted with a bold- 
face. We also define 



X 



12 



'12 



X 



1 ^21 
{x\ - 4)/6, 

{x' - x\)/rA, 



X12 

ni2 = n 
UA = riA 

where we follow the conventions of [2 



with 



X, 



(28) 
(29) 
(30) 

(31) 



We further assume these objects are on a quasi- 
circular orbit with leading-order angular velocity uj = 
{l/b){m/by^^ and orbital velocity V = {m/by/'\ The 
orbital separation b is assumed constant, as its time- 
evolution is driven by GW emission at high-order in v/c. 

The GR spacetime metric for such a binary is expanded 
as in Eq. (jl3p . In the near zone, the metric perturbation 
is given by 



^00 



2Ui + {1 ^ 2) + Oiv^) , 
2UiS,j + {1 ^ 2) + Oiv^) ^ 



(32) 
(33) 
(34) 



where 0{v^) stands for an (A/2)PN remainder, i.e. a 
term of 0{{v/c)^), and the notation -|-(1 o 2) means 
that one should add the same terms with the labels 1 
and 2 interchanged. The potentials Ua and Vai with 
A = (1,2) arc defined as 



Ua 



Pa 



VA^ 



(35) 



where pA and v\ = x\ are the density and the center of 
mass velocities of the respective objects, with the over- 
head dot standing for time differentiation. Field variables 
associated with a prime, e.g. p^, are to be evaluated at 
x' . In the point-particle limit, the metric becomes 



2mi ,^ 

/loo = ^ + 1 o 2 , 

ri 

Anil , 

ho^ -v\ + {1^2) 

n 

^^3 = '-y^^^3 +(1^2), 



(36) 
(37) 
(38) 



with remainders of relative O(v^). We have kept the PN 
leading terms in the metric that are proportional to niA 
only, but higher-order terms can be found in [s^ . while 
terms proportional to the spin of each BH can be found 
in [H. 



IV. SCALAR FIELD EVOLUTION 

In this section, we solve the evolution equation for the 
scalar field both for field points in the far and near-zones, 
as defined in Sec. IIV Al The former will allow us to eval- 
uate the energy flux carried by the scalar field at infinity, 
while the latter will be essential to find effective source 
terms that reproduce the known strong field solutions 
and to solve the evolution equations for the metric defor- 
mation. 



A. Zones 

As shown in Fig. [5J let us decompose the geometry into 
three regions: an inner zone (IZ), a near zone (NZ) and 
a far zone (FZ); see e.g. [SJ-jSg for further details. The 
IZs are centered at each object with radii TZiz- These 
radii are defined as the boundary inside which either 
r^'"J" 7^ or the usual PN approximation breaks down 
due to strong-gravity effects. We here take them to be 
sufficiently larger than niA and much less than b. The 
NZ is centered at the binary's center of mass with radius 
7?.NZ and excluding the IZs. This radius is defined as 
the boundary outside which time-derivatives cannot be 
assumed to be small compared with spatial derivatives 
due to the wave-like nature of the metric perturbation. 
We here take this boundary to be roughly equal to Aqw, 
where Agw denotes the GW wavelength. The FZ is also 
centered at the binary's center of mass, but it extends 
outside 7?.Nz- 

One can only apply the PN formalism when the grav- 
itational field is weak and velocities are small. When we 
deal with strong field sources like BHs and NSs, there- 
fore, one can use the PN scheme in the NZ and FZ only. 
In the IZs, one may not be able to use PN theory, since 
the gravitational field may be too strong. In this case, 
we have to asymptotically match our PN solution in the 
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FIG. 2: We consider three zones, inner zone (IZ), near zone 
(NZ) and far zone (FZ). The IZs are centered at each object 
and their radii TZiz satisfy TZiz <C h. The NZ is centered at the 
center of mass of the two bodies and the radius TZm satisfies 
7?.NZ ~ Agw, where Agw is the GW wavelength. 



can safely neglect the contribution from the FZ, since the 
fall-off of the source term is sufficiently fast. 

The solution in Eq. PO]) can be simplified by integrat- 
ing by parts several times and using that □[/ = —Anp 
and \3\x — x'\~^ = —At:5^^\x — x') to obtain 



d = I67r^ 



(fix 



M 



\X — X' 



+487r^ / P- 

P J M 



X 



X — x'\ 
1 

\x — x' 



M 



\X — X' 



d-'x' 



, as 



M 



U[,U[,5'^^\x-x')d''x' . 



(41) 



Expanding this solution in terms of particles 1 and 2, we 
arrive at 



(42) 



NZ with the strong field solutions valid in the IZs, inside 
some buffer regions that overlap both NZ and each IZ 
(see Refs. [Ill, [s^, for a description of how to carry 
this out in GR). The strong field solution for BHs was 
found in Refs. [28| and @ in the class of theories consid- 
ered here. 



B. Near zone solutions 



with 



, 167r ion/ nJPx' 



,a3 



M 



1 



\x — X' 



-4^C/i,,[/i,, + (1 2) 



d^x' 



(43) 



Since the NZ is in the weak field regime, we can apply 
the PN formalism to compact binary systems. Let us 
consider the even and odd-parity sectors separately. 



1. Even-parity sector 



The evolution equation for the even-parity sector is 
□,^ = - 647r2^p2- 647r2^p2 

P P 

- ^ (h^fi^i.M^^^^'' + h^fi^^.M''^''^^) , (39) 

with p = pi + p2 and remainders of Oih?'). 

First, let us consider weakly-gravitating objects, 
i.e. not BHs or NSs, in which case the PN expansion 
is valid also in the IZ. By substituting the GR PN met- 
ric of Eqs. the NZ solution to the above wave 
equation at leading PN order becomes 

g i« "1 /" /2 d^x' , 1^ "2 /■ ,2 <Px' 



M 



M 



d^x' 



\x — X 



^'(40) 



again with remainders oi 0{h^), with U = Ui + U2 and 
Ai denoting the constant-time, NZ+IZ hypersurface. We 



and 



03 

/3 



M 



{p'iUi, + p'2U[^,) 



\X — X' 



d'x 



(44) 



i?,„,f is the part of ■& that can be evaluated by considering 
a single object only, while t^croas is the part that depends 
on the fields of both bodies. 

The integrals that define both ■i^^.^if and z^croaa have sup- 
port in the IZs only, and thus, the NZ integral operator 
is homogeneous (source- free). When we discuss the NZ 
behavior of fields associated with compact objects, such 
as BHs or NSs, we cannot directly evaluate such IZ inte- 
grals. These are derived under the assumption that the 
PN expansion is valid everywhere, which fails for com- 
pact objects in the IZs. Instead, we need to determine 
these homogeneous solutions through asymptotic match- 
ing. Before doing so, it is helpful to study the meaning 
of each term for weakly-gravitating objects. 

Neglecting the size of the weakly-gravitating objects, 
the first term in Eq. P5|) in the NZ is evaluated as 



12 
Pi - 



d^a 



M 



(45) 
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with remainders of relative Oim/r), while the second 
term becomes 



0. 



' 1 JM 



M 



X' - y\ 



(46) 



The last equality can be shown by exchanging the inte- 
gration variables'^. Thus, one can approximate iJaeif as 



L'sclt 

ri 



4^ + (l-2), 



with the scalar monopole charge defined by 

P Jiz 



(47) 



(48) 



with A = (1,2). Here we put "IZ" to the integral to 
emphasize that the integration can be restricted to both 
IZs because the integrand is localized. 

The first term in Eq. (|47p represents the monopole field 
around object 1. These monopole fields give the leading 
PN contribution in the NZ unless both monopole charges 
qi and (72 vanish. This is indeed the case in EDGE the- 
ory, where (ai,a2,ck3) = (1, ^ 4, l)aEDGB- We will later 
show that this cancellation does really survive even if 
we consider NSs. If this cancellation occurs, the higher 
order terms of ©(m^/r^) in the expansion of Eq. (|45|) be- 
come the dominant contribution to ?9. The second term 
in Eq. (jlT]) is much higher PN order compared with the 
first term and hence sub-dominant in the NZ. 

Let us now consider 'd^^it for compact objects, where 
the IZ integrals must be treated carefully. Since the PN 
expansion is no longer valid in the IZ, one cannot use the 
simple extrapolation of the above result. In Sec. lIV Cl we 
match the NZ solution to the one obtained for isolated 
BHs in the strong-field 0, [2I] . We will not discuss the 
matching for NSs in this paper, but the order of magni- 
tude estimate 



the NS scalar monopole charge vanishes independently 
of the equation of state. Mathematically speaking, this 
is because the monopole charge is given by the integral 
of the Gauss-Bonnet invariant Rqb = — 4:Rf^^R^^'^ + 
RfiiypcrR^'''^'^ , which vanishes for any simply-connected, 
asymptotically fiat geometry. A more explicit proof is 
given in Appendix 1X1 

Let us now return to the t^crona contribution and 
consider first weakly-gravitating objects. To evaluate 
Eq. ([33]), one can use point-particle expansions of the 
potentials and the density, i.e. pA = niAS^^^ {x — xa) and 
Ua = TnA/TA- Simple substitution leads to 



Q;3TOim2 



n\ 



"-12 
62 ■ 



62 



O 



m 



(50) 



The first term in parentheses comes from the term 
Ui^iU2,i in Eq. (|44p . The remaining two terms come from 
the integral in Eq. (|331)- The second and third terms 
in parentheses look like scalar dipole moments for bod- 
ies 2 and 1 respectively. However, a Taylor expansion 
about the center of mass of each body, shows that the 
l/r\ piece of i^^rosB cancels, which implies that there is 
no scalar dipole. 

Let us now consider iJcroaa for compact objects. As 
discussed in the previous paragraph, one might expect a 
scalar dipole charge induced by the acceleration of object 
1 due to the gravitational field of object 2 (cx U2,i{xi)). 
In GR, however, acceleration is understood as geodesic 
motion in a perturbed geometry. The deviation of the 
local geometry from the unperturbed isolated geometry 
originates due to tides, and this is a relative 4PN effect. 
This is much smaller than the scalar monopole charge 
contribution from -d^^it. 

To summarize, the dominant contribution to 1? comes 
from the monopole charge associated with each object, 
which depends on its internal structure. 



E 



qi,A 



m , 



^ (49) 



should still be valid, where Ra is the radius of the Ath 
NS. When ai + a2 + Sa^ = 0, the cancellation observed 
in the weakly gravitating objects may still persist even 
for NSs. However, the cancellation will not in general be 
exact, except for the EDGE subcase. In EDGE theory. 



^ In fact, this integral vanishes to all orders in x. This is because 
is spherically symmetric, and thus, when it acts as a 
source to a wave equation, the solution should either scale as 1/r 
or it should vanish identically. We have here shown that there is 
no 1/r part. 



2. Odd-parity sector 

In the odd-parity case, the scalar field evolution equa- 
tion is 



(51) 



plus terms of 0{h^). Again, we first consider weakly 
gravitating objects. At leading PN order, the above equa- 
tion becomes 



Orjli) — -^-Cijk{hoo,mihkOJm + holjmhkl^im) 
— 32-^6jj/gt/^2m^^,jm J 



(52) 
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with remainders of relative 0{v'^). As in the even-parity 
case, we write the solution to this wave equation as 



(53) 



where 



P Jm 



1;^^ + (1^2), 



and 



(54) 



8^4 f ,w d^x' 



(55) 

Let us first consider self-interaction terms ^?aeit- Inte- 
grating by parts several times, we find 



1 



M 

U'i,p'iv'ik 



1 



M 



X — X'\ 

(fx' 



+ / Ulyi^J^^\x-x')d^x' + {1^2) 



M 



,(56) 



where we have used the relations DC/i = — 47rpi, \3Vi = 
-Airpiv'i and a\x - x'f'^ = -AttS'-^^x - x'). The third 
term vanishes when we take the point- particle limit^, 
i.e. PA = uiaS^^^x - xa), Ua = mA/rA and Vai = 
ruAVAi/rA- 

Let us evaluate the first and the second terms in the 
NZ. Keeping only the leading PN term in the NZ, we find 



(57) 



= 16^e.,fc^ / pU^fc., - Ul^v[,)d'x' + (1 o 2) 
P Jm 

= -f (1 O 2) , 

where we have defined 

,{A) 



(58) 



This leading-order PN term in -d^^^f represents a magnetic- 
type dipole. 

As in the even-parity case, to extend this result to com- 
pact objects we have to determine the value of pif^^ by 
matching the NZ solution in Eq. (|57|) to a strong field 
solution. This will be carried out in Sec. IIVCI for the 



BH case. For NSs, we just present an order of magni- 
tude estimate based on a simple extrapolation of weakly- 
gravitating results: 



(59) 



where S*^ is the spin angular momentum of the object. 
Following the procedure in Appendix \^ we can show 
that NSs cannot have scalar monopole charge in the dy- 
namical CS case. 

Next, we consider the cross term 'd^^^^ in the weakly- 
gravitating case. Integrating by parts several times, we 
find 



M 



I Ml 

2k.] 



1 



M 



d^x' 



d^x' 



+ / U'VL^5^^\x-x')d'x' + {1^2) 



M 



(60) 



One can take the point-particle limit of this expression 
without any trouble to obtain 



a4mim2 
-Id — - — ; — ei,fcUi2fc 



'■12"-l 



'-12"'2 



1 ^ n-^ 



' 1'2 



(61) 



Since plY^^ in Eq. (j57p do not depend on f , these terms 
are of relative 0{v^) compared to the leading-order term 
of^.„„. 

As for compact objects, the results found in the even- 
parity case also apply here. Terms proportional to 1 /r\ 
in the above expression suggest that each object has a 
dipole component induced by the companion. When we 
expand this expression around <g; 6, however, the 
terms proportional to l/(r^6^) cancel each other, as in 
the even-parity case, leading to no induced dipole mo- 
ment. Even if this were not the case, however, the cor- 
rections to the dipole moment would be higher order than 
the contributions from t^aeit- 

To conclude, the dominant contribution to § is clearly 
that of oif given in Eq. (|57p , which again depends on the 
structure of the source and thus violates the effacement 
principle. 



* The vanishing of this term is a general consequence of the sym- 
metry of the system. The source term contains an e'^* symbol, 
which must be contracted with other vectors to produce a scalar. 
We here have only two possible vectors to contract with, i.e. the 
velocity v\ and the unit vector n\ from object 1. Hence, any 
contraction with the Levi-Civita symbol should vanish. 



C. Matching near zone and strong- field solutions 
and finding the effective source terms 

In alternative theories of gravity, the point-particle 
limit is not always valid and the multipolc moments of 
compact objects may depend on the internal structure of 
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the source. In the previous subsections, we found that 
the dominant contributions to the scalar field come from 
self-interaction terms, which in turn depend on certain 
structure constants. In this subsection, we determine 
these constants by matching the -d solution to that of an 
isolated BH. 



1. Even-parity sector 

In the even-parity case, the monopole charges qi and 
q2 in Eq. (jTZ]) must be determined by matching to a BH 
solution. An isolated BH sources a scalar field [28[ , whose 
leading PN behavior is 



2q;3 mA 



Matching this solution to the NZ solution of Eq. 
obtain 



2a;3 

PmA 



(62) 



471) we 



(63) 



Notice that this monopole charge does not depend on 
(ai,a2), as for pure BH spacetimes, these coupling con- 
stants appear in combination with the Ricci scalar and 
tensor, which vanishes. This is to be contrasted with the 
NS case, in which qA depends on ai and a2 as well as 
and vanishes in EDGB theory. Interestingly, BHs do not 
have scalar hair in more traditional (Brans-Dicke type) 
scalar-tensor theories, while NSs do possess them. This 
situation is reversed in EDGB theory. 

The matching carried out above dealt with the 
monopole part of "d. That is, we have ignored any tidal 
deformation of either BH induced by its binary com- 
panion. In BH perturbation theory, one can calculate 
the deformation of the isolated BH metric to find that 
it depends on the sum of electric and magnetic tidal 
tensors, leading to a metric deformation that scales as 
{ri/hf {m2/h) for ri < 6 IMl IH-ill. Thus, in the 
IZ of object 1, tidal deformations lead to corrections of 
C(m^/6^), which are much smaller than the effects con- 
sidered here. Therefore, it suffices in this section to con- 
sider an isolated BH when matching the scalar fields. 

With this at hand, we can now treat BHs in even- 
parity, quadratic modified gravity as delta function 
sources of matter energy density, and with effective scalar 
charge density 



Pi) = qA 



X - Xa) 



(64) 



In the PN expansion such sources reproduce the BH so- 
lution found by Yunes and Stein at leading order. 

Let us make a few observations about the effective 
source term approach. First, notice that the scalar field 
diverges as rriA 0, which violates the small-coupling 
approximation. This is related to the fact that as one 
shrinks a BH, the radius of curvature at the horizon also 
goes to zero, probing increasingly shorter length scales. 



When the small-coupling approximation is violated, one 
can no longer neglect the scalar field's stress-energy ten- 
sor and the (^^^^ jT^ii/, /C^ji/) tensors that would dom- 
inate over the Einstein tensor. Of course, one cannot 
take this limit seriously, as we are considering here a low- 
energy effective theory, which is missing higher-curvature 
terms that would need to be included. Notice also that 
this is different from the behavior of scalar fields in tra- 
ditional scalar-tensor theories, where the scalar field van- 
ishes in the rriA — > limit. 



2. Odd-panty sector 

In the odd-parity case, the dipole charges of the respec- 
tive objects in Eq. ([57| are to be determined by matching 
against the appropriate BH solutions. An isolated non- 
spinning BH in the odd-parity case docs not support a 
scalar field. By contrast, a spinning BH does, and in the 
slow- rotation limit, neglecting higher order PN correc- 
tions, it is given by @ 



5 "4 i i 



(65) 



where XA = S^^/m^ is the normalized spin angular mo- 
mentum vector of the Ath BH. Matching this solution to 
the NZ ■d,,i! in Eq. ^7}, we obtain 



j _ 5 a4 j 



(66) 



With this at hand, we can now treat BHs in odd-parity, 
quadratic modified gravity as distributional sources of 
matter energy density and effective scalar charge density 



-p.\5^^\x-XA). 



(67) 



In the PN expansion, such sources reproduce the BH so- 
lution found by Yunes and Pretorius [6[ at leading order. 

Let us make a few observations about this solution. 
First, notice that the pseudo-scalar dipole charge is well 
behaved in the limit ruA — >■ 0, because there is a max- 
imum BH spin Ix^l < 1. Second, notice that in the 
\xW ~^ limit, this dipole charge vanishes, which is a 
consequence of Birkhoff's theorem holding in CS grav- 
ity [H, 13 113 ■ Namely, non-spinning BHs in CS the- 
ory are the same as BHs in OR (i.e. Schwarzschild 
BHs). Therefore, in this case the point-particle limit is 
well-justified and the metric deformation or the scalar 
field does not depend on the internal structures of non- 
spinning sources. 



D. Far-zone field point solutions 

Let us assume that wc have the wave equation 

Ur,d = T{t,x), (68) 
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where r denotes the source term. The far-zone field point 
solution to this wave equation is given as l63j 



^ oo 

-T 



47r ^ — ' m\ 

7n—0 



-dm 



M 



T(U,X )X 



, (69) 



with u = t — r. By using u^i = —rii and by keeping only 
terms proportional to 1/r, the above solution reduces to 



11^19" 



E 



An r ^ ml dV' 

m=0 



M 



(70) 

Here, the region M denotes the hypersurface of t — r = 
const. In the following, we apply these formulas to the 
even and odd parity cases separately. 



1. Even-parity sector 



Following the discussion in Sec. IIVB| the evolution 
equation for the scalar field is dominantly 

□^7? = -47rgi(5(3)(a;- cci) + (1 2). (71) 

From Eq. (|70p . this wave equation can be solved as 



1 d'' 



r ^-^ m! dV ' i m 
+ (1 o 2). 

The TO = term gives 



S^'^Hx' - xi){n, x'^rd^x' 



(72) 



(73) 



where we have defined the total scalar monopole charge 
q = qi + q2. Recall that this monopole charge q refers 
to the scalar field, and not to an electromagnetic one. 
For a BH binary or a NS binary in a quasi circular orbit, 
q only changes during merger, as mass is carried away 
in radiation. Thus, monopole radiation is inefficient and 
suppressed. 

For the m = 1 case, we find 



(74) 



where we have defined the total scalar dipole moment as 
D' = qix{ + q2x\ . (75) 
When we evaluate this for circular orbits, we find 

^- = iU^-..^).i.:nN (76) 
r \ m TO / 

where we have defined the relative velocity = Wi — V2 ■ 
The TO = 1 term clearly leads to dipole radiation in the 
FZ, which is less relativistic than GR quadrupolc radia- 
tion, becoming stronger at smaller velocities. Of course. 



this term is proportional to the coupling constants of 
the theory, which are assumed much smaller than one. . 
Reference [l^ failed to recognize such dipolar emission 
because they considered the motion of test particles that 
had no scalar charge. We cannot think of any mechanism 
that would suppress such dipolar radiation. 



2. Odd-parity sector: spinning bodies 

As in the previous Section, the evolution equation for 
the scalar field is dominantly 



= A^^^l\^'•^\x - Xi),, + (1 o 2) , 



(77) 



By using Eq. (|69p . the far-zone field point solution is 
obtained as 



^ — ^ m ' 



m=0 



fi\S^^\x' -xi).,x'^^d'^x' 



M 



+ (1 o 2) 



(78) 



When TO = there is obviously no contribution to the 
scalar field. When to = 1, 



[ S^^\x-xi)^,x'd^x 
Jm 



and thus 



7?" 



(79) 



(80) 



with fii = /iii -f fi2i- Notice that we recover the solution 
of Yunes and Pretorius 0] for the first term of the above 
equation with given as in Eq. (j66p . These terms will 
not strongly radiate because /i^ is non- vanishing only for 
spin-precessing systems. Even then, such radiation would 
be suppressed by the ratio of the orbital timescale to the 
precession timescale. 

The TO = 2 contribution, by contrast, depends on the 
much shorter orbital timescale. We look for terms of 
0{r^^) since they are the only ones that contribute to 
the energy flux at infinity. Keeping in mind that the 
function being differentiated depends on retarded time, 
we can rewrite Eq. ((78)) as 



i9 



M 



r ' — ' to! dt 

m 

+ (1^2). 
When m ~2, we have that 



^1,5(3) (a;'-a;i) ,;(nfc x"')"'d^x' 



(81) 



Ail / 5^^^ (x - Xi),,xPxid^x + (1 o 2) = -2^iP'^ , (82) 
Jm 

where the pseudo-tensor quadrupole moment (not to be 
confused with /^ V"*) is defined as 



{i j) . {i j) 
^1 f^l "1" ^2 



(83) 
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The m = 2 contribution becomes 



r 



(84) 



where the final equality is evaluated on a circular orbit. 
Notice that such a scalar field will strongly radiate be- 
cause ^^-^ depends on the orbital tiniescale. 



case, there cannot be any FZ vectors n' present. Thus, 
all vectors that can be contracted onto the Levi-Civita 
symbol must lie in the same orbital plane and this obvi- 
ously vanishes. This argument should be true at all PN 
orders^. 

Let us then consider the next-order term. This will 
arise from the leading-order source term [right-hand side 
of Eq. with m = 1 in the NZ sum: 



3. Odd-parity sector: non-spinning bodies 

When both objects are non-spinning, the self- 
interaction terms produced by the effective source iden- 
tically vanish. One is then left with the source term 
constructed from the product of the gravitational fields 
of objects 1 and 2. These terms will be proportional to 
miTO2- As we will see, there are many contributions that 
turn out to vanish upon NZ integration. For pedagogi- 
cal reasons, we will show here explicitly how this happens 
and eventually arrive at contributions that do not vanish. 

The evolution equation for the scalar field to leading 
PN order is 



a4 



32—eijk'mim2Vi2k — — 



1 



1 



(85) 



where we substituted the NZ metric components in the 
point-particle approximation. The leading order term of 
the solution to this differential equation, i.e. the m = 
term in the sum of Eq. ([70]) , is evaluated as 



8 a4 



t^Fz = -^rnim2e,.jk— / ( — ) ( — 



1 



= „i6^^,^,e,,,^af )af )C^9(?)y ^ o .(86) 

p J, J 

Here we integrated over the NZ-I-IZ hypersurface 
without taking any care of the strong gravity region in 
the IZs. One can easily show that the contribution from 
the IZs is not large in the present case. In the second 
line, we replaced partial derivatives with respect to 
acting on 1 /rA with (minus the) particle derivatives with 
respect to x^: 



d 



d 



-d: 



(A) 



(87) 



with A = (1, 2). We commuted these particle derivatives 
with the integral, and finally obtained a typical NZ inte- 
gral, discussed in Appendix IbI From Eq. (jB4p . we know 
that Y — h, and by taking all particle derivatives, the 
last equality is established. 

We could have inferred that the m = term in the 
sum does not contribute for non-spinning BHs without 
any explicit calculations. The argument here is similar 
to that in footnote H] Possible vectors to contract with 
the Levi-Civita symbol include the velocities v\ and the 
unit vectors n\ i but not spin vectors S\ , as we here con- 
sider non-spinning BHs. In particular, for the m = 



_ 8 q;4 mim2 p ^ 

^FZ — _ ^ ^ ^ ^ijh^l2 
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d 



dt Jm V '"1 



1 



= -16 



a4 TO 1 7712 



r 



dt 



where we have used Eq. (|B3p . which defines Yp. By direct 
evaluation, one can show that this term also identically 
vanishes. The first non- vanishing contribution coming 
from an TO = 1 term must then be 0{v^) smaller than 
the ordering of the to = term. 

Finally, let us consider the (next)2-order term. This 
can arise only from the leading-order source term with 
TO = 2 in the NZ sum: 



4 a4 TO1TO2 
TT f3 r 



^12 



X 7'2 , 



Q!4 miTO2 



„P9f .., -9(1)5(2)5(1)5(2) 



Y, 



which simplifies to 



ai 1 r]mSm 

1 — 

pro 



knipU!'^v^2'rii2 ' 



(89) 



(90) 



where we have defined the mass difference Sm = mi — TO2 
and the symmetric mass ratio r] = mim2/m?. We have 
here used Kepler's law and expanded the STF tensors. 
This is the dominant FZ behavior of the scalar field, 
which as we see is much suppressed relative to the odd- 
parity solution we found for spinning BHs. 



E. Summary of this section 

Let us summarize the results found so far for later use. 
In the even-parity case, generically at least one of the 



^ One may think that one can construct a vector that does not he 
in the orbital plane by taking the cross product of two vectors 
that lie on this plane, e.g. n\2 X v\2- However, since GR is 
parity even, such a vector cannot be present in the PN metric. 
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A 


B 


c 


a 


6 c 


d 


Even-P 


qi 







1 


- 




Odd-P, Spins 


i i 





ii 


2 


- 




Odd-P, No Spins 








2 


2 2 


2 



TABLE I: Scalar field parameters, as defined in Eqs. (|9H) 
and (|93p . The quantities q\ and are defined in Eqs. (|63|) 
and (|66p . while aj^^ is defined in Eq. (|92p . The quantities Di 
and /^i are defined in Eqs. H75|) and H78|) . while (jpl is given 
in Eq. (fMj) . 



binary component objects will have a scalar monopole 
charge. Since the scalar field excitation due to the in- 
duced monopole is dominant, we neglect all the other 
less important contributions. Weakly gravitating objects 
need not have a scalar monopole charge if ai+a2 + ici^ = 
0; and BHs have no scalar monopole charge if as = 0. 
In EDGB theory, NSs have no scalar monopole charge. 
In the odd-parity case, the dominant contribution is the 
magnetic-type scalar dipole moment induced by spins. 
Generically, astrophysical objects will possess spin, but 
we will continue to include non-spinning results to com- 
pare with previous work. 

In the NZ, we can parametrize the leading PN terms 
of the scalar field as 

^- = 4^ + l|l + (1^2), (91) 

where (A, i?, a, 6, c, d) are given in Table U and for com- 
pactness of the Tabic we define 

all ^ -lQjW^v,sv{2 ■ (92) 

In the FZ, we can parametrize the scalar field as 

^9pz = - , (93) 
r 

where C is also given in Table U and we define 

2 

all ^ 16jVmSm'^e,,,v',^n{l , (94) 



V. METRIC EVOLUTION 

In this section, we solve the evolution equations for 
the metric deformation in the FZ, so that we can calcu- 
late the gravitational energy flux at infinity. Note that 
throughout, we use the Newtonian relationship = in/b 
(and similarly for the acceleration). This relationship 
must be corrected at higher PN order or at ©(<;). As we 
mentioned earlier, here we do not take into account the 
corrections to the orbital motion due to the conservative 
force at C(^). These conservative effects do not interfere 
at ©(t^) with the radiative effects that we are concerned 



with in this paper. Therefore the corrections to the GW 
waveform become a simple summation of these two dif- 
ferent types of effects. 

For the FZ field points, the solution to the metric de- 
formation equation of motion [Eq. (jl7|) ] can be read from 
Eq. dZni): 

= --rT.^^^ 4 in^."'rd^.' + O (.-) , 

(95) 

where we have defined the source term as 

^a,{^J^+jS^+a,K^-\r^K (96) 

Notice that this corresponds to an IZ-I-NZ integration for 
FZ field points, where we have neglected the FZ integra- 
tion because it is subdominant. 

The integrals presented above have to be carried out 
also in the IZ, where the PN expansion is not valid any- 
more. In GR, however, such divergences can be ignored, 
using a regularization scheme. Since both the true so- 
lution and an appropriately regularized solution satisfy 
the field equations in the NZ, their difference due to the 
IZ contribution is only through a homogeneous solution. 
Such homogeneous solutions are regular in the NZ and 
FZ, but can be divergent in the IZ. They are character- 
ized by the multipole moments of the respective objects, 
which can be determined by studying tidal perturbations 
around a strongly gravitating object. One can then per- 
form matching of the metric solution, as for the scalar 
solution, but the metric matching is beyond the scope of 
this paper. In what follows, we only consider the regular- 
ized contribution, following Hadamard partie finie (FP) 
regularization |64{ . We comment more on the divergent 
contribution at the end of this Section. 



A. Even-parity sector 

Let us focus on the metric perturbation in the even- 
parity sector first. The leading order term both in the 
PN and 1 /r expansion at infinity is formally given by 

O^.^l j^Tlpd'x, (98) 
f,^^-^ / J,,d'x, (99) 

where we have defined J,j = ^jj"^ + J^p . The source 

terms H^i^ and do not contribute to this expression 
since they identically vanish in the NZ where R^i, = 0. 
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We can estimate the order of magnitude of both t) 



J 



and f)^ as follows: 



J 



r ^ 



(100) 
(101) 



Here we factored out in the final expressions, since the 
GR leading quadrupolar field is also proportional to . 
Clearly, the dominant contribution comes from Eq. pOOp . 

Let us now make this computation more precise. The 
stress-energy tensor will contain self-interactions of the 
form '&A.i'&A,j and cross terms of the form 'i?i,i??2j"- The 
former case leads to divergent integrals, which must be 
determined by strong-field matching, so we do not con- 
sider them here. Let us concentrate on the latter, which 
take the form 



1 



(102) 



3) 



(103) 



which sources the metric perturbation 
4 



^ J M 

An 



^(l)j)(2) 



(lo2), 
(104) 



where we used an integration formula for the triangle 
potential given in Appendix |B] We can see that this cor- 
rection is OPN relative to the radiative metric perturba- 
tion in GR, just as we predicted in Eq. (|100|) . However, 
this correction turns out to be still smaller in the energy 
fiux than the dipole scalar radiation, which gives a -IPN 
correction. 



B. Odd-parity sector 

We now focus on the odd-parity sector, for which 
the solution is given by the term proportional to a4 in 
Eq. namely 



Ujj — 



8q;4 



r 



M 



(105) 
(106) 



The stress-energy contribution f)^ is the same as in 
Eq. (fTUa . 



The K. contribution to Eq. (|105p is more involved. The 
leading-order behavior of the K, tensor is 







ira.lk "T fl'tk^irn J 
+^,fc;^j'ni(^iO,mfe + hrnk,iO — hkO,im ~ ^im.feo) 



^k^ikl{2hQ[mj]i,n — 2hi[j^„i]jn ■ 

'2'deikih^,y^,n]i,n + {i ^ j) ■ 



hoo,ji ) 



(107) 



Other terms are of higher PN order. By applying the 
Lorenz or harmonic gauge condition h^'^ .v = 0, substi- 
tuting hij — hooSij into Eq. (|107|1 . and using ejkih^'^'''^ = 
f^jki-^M 0, we get 

IC[^^ — 2d ^k^ jklhoo.il ~ '^'^ ,km^jklho[rn,i]l 

— "^"^ .kf-]klho[m.i]lm + .k^jklhoO,il + (i j) ■ 

(108) 



The ]Cij term in Eq. (|105p is then a sum of four terms, 
namely 



(109) 



where we have defined 





16a4 
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Jm 




^ 16a4 
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JM 




^ 16a4 
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Jm 


W = 


16a4 








f 

JM 



'&,kejkihoo,ticPx + (i o j) , (110) 

■&,kTn<^jklho[m,i]ld^X -t- (i O j) , (111) 

^,kejklho[ma]lmd^X + (j o j) , (112) 

T^Mjkihao.ad^x + (i o j) . (113) 



When we substitute the PN metric into the above 
terms, the right-hand sides depend on the velocity vectors 
v\ (which depend on time only). The field ?9 is given in 
Eq. (|9ip and its derivative can be computed simply from 
that equation. Since this field is a NZ one, it depends on 
time through the positions of the objects, which implies 
that its time derivative can be converted into a spatial 
derivative via dtf{ri) = —v\dif{ri). 

Let us begin by making a simple order of magnitude 
estimate of how large the regularized contribution is. For 
this, it suffices to look at Ec^s. pITO)) and PTU|) : 



\m r 



The 1? field here is that of the NZ, and hence 



r ^ 



(114) 
(115) 

(116) 
(117) 
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where x stands for the magnitude of Xi a-nd X2- From this 
analysis, f)^ is clearly larger for rapidly spinning objects, 
leading to a 2PN effect. 

For the non-spinning case, one might expect the K, 
contribution to lead to a 6PN effect, but as we explain 
in Appendix [Cl these leading-order effects actually van- 
ish. This cancellation can also rather easily be seen by 
integrating by parts in Eqs. (jll0p - (|113|) . After discard- 
ing boundary terms (taking into account the boundary 
term is equivalent to adding homogeneous solutions, cor- 
responding to deformed multipole moments of compact 
objects), we obtain expressions of the form ejki "& /ioo,fcz...j 
which obviously vanishes by the antisymmetry of the 
Levi-Civita tensor. We carry out a more careful analysis 
in Appendix[Cl where we explicitly show that the leading 
and first sub- leading order terms vanish^. The first non- 
vanishing term is then of 0{v'^) smaller than the order 
of magnitude estimates in Eqs. (|116p and (|117p . leading 
to 7PN and 4.5PN contributions at 0{x^) and O(x^), 
respectively. 

Since the largest contribution seems to arise for spin- 
ning BHs from the t)^ term, let us consider this in more 
detail. Two possible contributions are generated here: 
one that depends only on self-interaction terms, and one 
that depends on the cross-interaction. The former leads 
to divergent integrals, which need to be matched from 
strong-field solutions, and we do not consider these here. 
The latter leads to the metric deformation 

IS - - v^^'^^ {^d[]}dfY - + (1 o 2) 

+ 3<2 [5 ("^2MuO ("'12/^2;) - f^lkf^2] } + (1 ^ 2) , 

(118) 

which is clearly of the order predicted in Eq. ()116p . 
i.e. 2PN order relative to GR. This is of the same order 
as the energy flux correction carried by the pseudo-scalar 
radiation. 



C. Multipole moments 

In this Subsection, we discuss the additional contribu- 
tion from the IZs, which enter as additional homogeneous 
solutions in the NZ and FZ, These contributions are ho- 
mogeneous in the sense that they arise from sources that 
have support only in the IZs, and thus they vanish in 
the NZ and FZ (see e.g. the discussion prior to Eq. ps]) ). 
The homogeneous solutions are characterized by the mass 
and current multipole moments of the strong-field bodies, 
which must be determined by matching to strong gravity 



In Appendix [C] wc only show this for non-spinning BHs, but a 
similar calculation can be performed for spinning BHs to 0{x)- 



solutions in the IZ. When we solve the non-linear equa- 
tions of motion iteratively, the source terms in general 
can be classified into two pieces: a self-interaction part 
and a cross-interaction part, as in the case of i? in Sec. IIVI 
The cross-interaction part is sourced by the companion, 
while the self- interaction part is not. 

The self-interaction part is rather easy to handle be- 
cause matching involves only a single isolated object. As 
described in Sec. lIVI3l these self-interaction terms can be 
thought of as homogeneous solutions that have support 
only in the IZ. As such, in the small-coupling approxima- 
tion, they satisfy homogenous field equations that take 
Einstein form. If the spin of the object is neglected, the 
only possible linear perturbation to such a homogeneous 
solution that is compatible with asymptotic flatness is 
a shift of the body's mass (in the 1/r piece of the {t,t) 
and diagonal parts of the metric). In essence, this is a 
consequence of Birkhoff's theorem, which holds for ho- 
mogeneous solutions. Such a shift is consistent with the 
strong-fleld, non-spinning BH solution in EDGB theory 
found in [28|. In that case, the mass shift is simply 
niA ^ (49/80)C3TOA. 

For spinning objects, one expects there to be higher 
multipole moments in the strong-field solution. However, 
one should be able to absorb current dipolc moment mod- 
ifications by a redefinition of the spin parameter, while 
the mass dipole moment will be absorbed by the redef- 
inition of the position of the center of mass. Therefore, 
the leading-order corrections that survive are the mass 
quadrupole moment, which produces a metric perturba- 
tion in the NZ proportional to As we will see, when 
we consider FZ solution, there is an additional factor of 

that enters. 

Therefore, contributions to the energy flux from the 
quadrupole or higher multipole moments are at least 3PN 
order relative to that from the GR quadrupole formula. 
We will later find that corrections to the energy flux due 
to scalar radiation appear at -IPN and 2PN relative order 
for the even and odd-parity cases, respectively. Hence, 
the contributions from the multipole moments that we 
discussed here are deflnitely smaller than those intro- 
duced by scalar radiation in the even-parity case, and 
at most, the same order in the odd-parity case. 

Let us take a look at spinning BHs in the odd-parity 
sector in more detail. At 0{x) there is freedom in adding 
a homogeneous solution proportional to 1/r^ in the hoi 
component. This corresponds to a freedom in shifting 
the Kerr parameter measured at inflnity. Reference [6| 
set this homogeneous solution to zero so that there is 
no shift in the Kerr parameter. At 0{x'^), there should 
be corrections proportional to in hij which shifts 

the quadrupole moment. Since there is no parameter 
in the Kerr geometry that can absorb this correction in 
the quadrupole moment, this correction cannot be 
eliminated. 

The effective source term that reproduces this correc- 
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tion should look like 

□ f]y = -'^'KQiUiUj{5ki~'iSi,kSi,i)5^^\x-xx)M 

+ (lo2), (119) 

where Qa — 0{C,imAa\) and SA,k = S\/m\ is a unit 
spin angular momentum vector. The solution of this wave 
equation at 0{l/r) is given by 

1 1 f)™ 

m=0 

Jm 

(120) 

The leading-order contributions at m = (2PN) and 
TO = 1 (2.5PN) vanish, leading to the first non-zero con- 
tribution at TO = 2 

= O (^Qc^v) = a^v' X O (x'v^) , (121) 

which is 3PN relative to GR. Therefore, the self- 
interacting correction in the metric at 0{x^) is smaller 
compared to the corrections in the energy flux carried 
by the scalar field and the metric field with regularized 
modification. 

The cross-interaction part is more complicated. In this 
case, we have to consider the induced multipolc moments 
due to the presence of the secondary object. Thus, even 
if we consider non-spinning objects, higher multipole mo- 
ments might be induced. Another important difference is 
that neither the mass monopole nor the spin dipole can 
be simply absorbed by a redefinition of the mass and spin 
of each object. This is because the shifts of these mul- 
tipole moments depend on the orbital parameters, such 
as separation b. Notice, however, that the effects of the 
secondary object propagate only through the scalar field 
or the gravitational tidal force. 

The order of magnitude of the former scalar field ef- 
fect is more complicated to estimate and it depends on 
the situation. In the even-parity case, sourccd by the 
secondary body at the position of the primary body is 
proportional to 1/6. In EDGB theory, since i? has shift 
symmetry within the context of the classical theory, the 
effects are suppressed by the gradient of the field, i.e. they 
are proportional to l/&^. In the odd-parity case, there 
is again shift symmetry and the monopole scalar charge 
is absent. Because of these two reasons, the suppression 
is proportional to in CS theory. These suppressions 
will be sufficient to conclude that the effects are relatively 
at least IPN and 3PN in the even and odd-parity cases 
respectively, which is smaller than the effects induced by 
scalar radiation. 

In the odd-parity non-spinning case, the latter gravi- 
tational tidal force dominates over the scalar propaga- 
tion effect. To calculate this tidal force properly re- 
quires asymptotic matching between the IZ solution and 



a strong-field, perturbed Schwarzschild solution in CS 
gravity. Perturbations of the Schwarzschild spacetime 
can be decomposed as a sum over electric and mag- 
netic tidal tensors (see e.g. [S^l). The former scale 
as 1/6'^(1 + V + v'^ + ...), while the latter scales as 
v/b^{l + V + v'^ + . . .) [5^. Such tidal deformations will 
induce gravitational waves that will scale as the second- 
time derivatives of the electric and magnetic quadrupole 
deformations, i.e. they will scale as uj'^ /b'^{l + v + v'^ + . . .) 
and uj^v/b^{l + v + + . . .). In GR, the leading order 
effect is induced by the electric quadrupole moment and 
it scales as uj'^/b'^, a 5PN order effect. In CS, we expect 
the magnetic quadrupole moment to provide the leading- 
order deformation, and at the level of energy flux, this 
couples to the GR metric perturbation produced by ci- 
ther the mass octupole or current dipole moment, leading 
to 6PN correction. This interpretation seems to be con- 
sistent with the results of Pani et al. . 

VI. ENERGY FLUX 

The inspiral of a compact binary system is controlled 
by the system's change in binding energy and angular 
momentum. The binding energy changes according to 
the dissipation of energy carried by all dynamical fields, 
which here includes the metric perturbation and the 
scalar field. The stress-energy tensor (SET) associated 
with each field quantifies the density and fiux of energy 
and momentum. The energy loss is calculated as the in- 
tegral of the energy fiux through a 2-sphere of radius r 
in the limit r 00 and in the direction of the sphere's 
outward unit normal n'. That is, for some field ip (be it 
/ly , , or z?) with SET tI£\ 

= lim / Ir'^fn') r^dn, (122) 

r-i-oo Jg2 \ / LU 

where the angle brackets with subscript lu stand for orbit 
averaging. 

The total energy flux can be ordered in powers of <r as 

E^Eoi, + ^SE + 0{<;^). (123) 

The GR energy flux E"^^ is given by the GR metric per- 
turbation only, without any contributions from the scalar 
field at 0{'^^), as there is no scalar field in GR. For cir- 
cular orbits, this is 

i;;,, = -^^V°. (124) 

The C(<r) correction, SE, can be decomposed into 

SE = SE^'^^ + SE*-^^ , (125) 

where the first term is the scalar field contribution and 
the second term is the contribution of the deformed met- 
ric perturbation. 



18 



The scalar field contribution is calculated with the SET 
given by Eq. ([7]): 



which leads to 



(126) 



r — > oo 



limit, I? must be that 



Since we are taking the 
vaHd in the FZ. 

The metric deformation contribution to the energy flux 
is slightly more subtle. This modification to the GR flux 
can have three distinct sources: (i) the effective SET in 
terms of fej, and f)ij may be functionally different, but as 
shown in [401, ^^i^ ^"^^ for the class of theories we 
consider here^; (ii) The orbital equations of motion, and 



the associated relations m/h 



and uj = V /m, 



be modified at C'(<^), as was partially calculated in j28l |: 
(iii) The generation mechanism of the FZ metric pertur- 
bation is modified, i.e. the radiative part of the metric 
perturbation is deformed. We consider here only the dis- 
sipative modifications introduced by (iii), as (ii) would 
require an analysis of the equations of motion, which is 
beyond the scope of this paper^. 

Letting Hap = hna + ^^ap + C'('^^), the effective SET 
of GWs is given by |40j 



1 

32^ 



rjTT LTa/S 



(127) 



where the angle brackets with a subscript A stand for 
a quasi-local average over several wavelengths and TT 
stands for the transverse-traceless projection 



rjTT 



A 



1 



ij.kl 



-P^jPki, (128) 



with Pij = Sij 



j the projector onto the plane perpen- 
dicular to the line from the source to a FZ field point. 
Expanding this SET in orders of the 0{^^) part leads 
to Egvi, while the 0[<;) part is 



167r 



(129) 



^ Reference [40l | showed that the TT gauge exists in quadratic grav- 
ity as r — >■ oo. Any non-TT propagating mode that is sourced 
in the NZ vanishes in the FZ at all orders. This is in contrast 
to scalar-tensor theories in the Jordan frame, where the scalar 
"breathing" mode is present in the metric. This difference comes 
from the way the metric deformation and the scalar field couple 
in the field equations. In the quadratic gravity case, i9 does not 
multiply Gfii/ in the field equations (the Einstein-Hilbert sector 
of the action is unmodified), while the opposite is true in scalar- 
tensor theories in the Jordan frame. Therefore, in the former 
f)/ji/ and i9 decouple in the r — > oo limit and there is no breath- 
ing mode. In contrast, in the latter the coupling between fi^,!/ 
and i9 remains in the limit r — > oo, leading to a non- vanishing 
breathing mode and a modification to the effective SET. 

* The distinction between (ii) and (iii) can be ambiguous at higher 
PN order, because how the orbital parameters arc modified de- 
pends on the gauge choice. However, as long as we impose the 
harmonic gauge condition on both GR and the deformed metric 
perturbations, we do not have to worry about this gauge issue 
at least up to next-to-leading PN order. 



(130) 



As before, the hag and Ijag are those valid in the FZ. 



A. Scalar field correction to the energy fiux 

1. Even-parity sector 

In the even-parity case, i!}^^ is dominated by the dipole 
component [Eq. (|74p ]. which we repeat here for conve- 
nience: "d^^ = Diu'^/r, where Di is the NZ dipole given 
in Eq. ([75]) . This is inserted into the energy loss formula, 
Eq. (jl26p . Since the FZ scalar field depends on retarded 
time, both time and spatial derivatives can be written as 
time derivatives of the NZ moments. This gives 

5E^'^^=-[3j (d.Dju'^) dn = -^(3(d'D,) , 

(131) 

which for circular orbits gives 

6 6 mr 

(132) 

Note that here, as before, the m — limit diverges, 
because the effective theory breaks down on short length 
scales and <j ^ 1 is violated. 

When the compact bodies are BHs, their scalar 
monopole charges are given by Eq. ((63|). qA = 
2a^/ {(imA), which then leads to 



3 if' rri^ 

Comparing this with the GR energy flux, we find 



(133) 



5E^ 

EGR 



5 1 5m^ 

T^Cs — — 
9o ?7* 



(134) 



a relative -IPN effect. That is, the energy lost to the 
scalar field due to dipole radiation would enter as a lower- 
order in V effect than the energy loss in GR. If one takes 
the limit m2 — >■ oo while keeping (mi,?;) fixed, then the 
above ratio scales as mj~^; i.e. the energy flux ratio is 
sensitive to the smallest horizon scale of the system. The 
effect is of a similar size for comparable stellar-mass bi- 
nary and EMRI system. A SMBH-SMBH binary experi- 
ences the smallest effect. 



2. Odd-parity sector: spinning bodies 

The scalar field 1}^^ is here dominated by the 
quadrupole component [Eq. (I84[) ]. which we repeat here 
for convenience = flijn^^ /r = —ui'^fiijn^^/r, where 



the quadrupole tensor fj,ij is defined in Eq. (1831) . Insert- 
ing this into the energy loss formula [Eq. (|126p ] gives 



J 



47r 
15 



/3 



(135) 



Let us evaluate this for quasi-circular orbits with 
non-precessing spins. The third time derivative of the 
quadrupole tensor /i.^ becomes 



and the total energy flux is 



(136) 



(137) 



where vi2 is the unit vector in the direction of the relative 
velocity and the dimensionless quantity A is defined as 



^ = Xl'^i X2J2 ■ 

m m 



(138) 



Notice that 5E^^'^ in Eq. is finite in the EMRI limit. 

Note also that when both spins are perpendicular to the 
orbital plane, A is as well, and the second term of SE^"^^ 
vanishes. Comparing Eq. ()137p with GR, 



5E^ 
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C,-^[A' + 2{{A-v,2r)Jv\ (139) 
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hence scalar radiation in the odd-parity sector is clearly 
a relative 2PN effect. This effect was not included in the 
work of Pani et al. [65|, who found a 7PN correction, 
since their simulations did not include spins. If one takes 
the limit m2 — > oo while keeping (jni,v) fixed, then the 
above ratio scales as mj~^m^^; i.e. the energy fiux ratio is 
sensitive to the geometric mean of the two horizon scales 
in the system. This implies that the effect is greatest for 
comparable stellar-mass binaries. 



3. Odd-parity sector: non- spinning bodies 




FIG. 3: Comparison of Eq. (|141|l to the numerical results of 
Pani et al. [63|. The latter can be mapped to the generic 
quadratic gravity action of Eq. Q by letting Q4 = — acs/4, 
which then implies that ^4 = — Ccs/16. We here used (4, = 
6.25 X IQ-'^, which is equivalent to their parameter (^cs ~ 0.01. 
Observe that at low velocities, in the regime where the PN 
approximation is valid, the two curves agree. 



which shows that this is a relative 7PN effect. In contrast 
with the cases of even-parity and odd-parity with spins, 
this effect is dominantly controlled by the total mass, 
rather than the mass ratio. The effect is greatest for a 
system of stellar-mass BHs. 

The above result can be compared to numerical cal- 
culations recently performed by Pani et al. 



65|. They 



estimated the effect of scalar radiation in dynamical CS 
gravity (ssj for non-spinning, circular EMRIs. They nu- 
merically solved the master perturbation equations on a 
Schwarzschild background to obtain the time evolution 
of the scalar field and the metric perturbation, caused by 
a non-spinning point particle. Figure [3] compares their 
results to ours, found in Eq. (|14ip . Observe that the nu- 
merical results of Pani et al. are in excellent agreement 
with our post-Newtonian calculation, which extends it to 
comparable mass-ratios (notice the factor of Sm/m). 



The odd-parity t9pz in Eq. ((90|) can be used to evaluate 
the energy loss in Eq. (|126|) : 



Se'^^) = _ 256KC4<5m^77^ ( ^ j J dfl 



256 

64 
15 



dt(^e'^''n'Pv'l2n{l 



Compared to the GW radiation in GR [Eq. (|124p ]. this 
scalar radiation becomes 



SM'^^ _ 2 Sm"^ 

iJGR 3 77l2 



Metric deformation correction to the energy 
flux 



For the even-parity case, the correction to the energy 
flux that arises from the deformation to the gravitational 
metric perturbation is at least of OPN order relative to 
GR. This is higher PN order compared to the scalar 
dipole radiation found in Sec. IVI A 11 and thus, we will 
not consider it further. 

For the odd-parity case with spinning BHs, one of the 
leading contribution comes from the metric correction 
sourced by T^P , which is given in Eq. (|118p . Inserting 

(141) 

this metric perturbation into Eq. (|130p . the energy fiux 
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correction relative to GR becomes 



Sm 75 C4 /a^aU0^ 



(142) 



which is of relative 2PN order, just as the contribution 
due to scalar radiation in Eq. (|139p . Notice that both 
the metric deformation and scalar field corrections to the 
energy flux are of O(x^), but the latter is larger by a 
factor of 0(77-1). 

We expect 0{x) corrections to the energy flux due to 
the metric deformation to be higher PN order. For very 
slowly spinning binaries, however, they may give larger 
corrections compared to the 0{x^) 2PN ones presented 
here. 

In the odd-parity sector with non-spinning objects, the 
regularized contributions to the metric deformation can 
only provide energy flux corrections of at least 7PN order. 
However, as explained in Sec. IV Cl we expect that match- 
ing strong-field solutions to the non-regular NZ ones may 
generate 6PN corrections in the energy flux, similar to 
those found by Pani et al. [65j . 



VII. IMPACT ON GRAVITATIONAL WAVE 
PHASE 

How do all these modifications to the energy flux af- 
fect the GW observable? To answer this question, we 
compute the Fourier transform of the phase of the GW 
response function in the stationary phase approximation 
(SPA), where we assume the GW phase changes much 
more rapidly than the GW amplitude (66j . 

We begin by parameterizing all the corrections to the 
energy flux that we have studied so far via the following 
power law: 



E = Eg^{1 + Av'') , 



(143) 



where {A, a) are summarized in Table HIl for the four dif- 
ferent sectors considered. 

With the generic energy flux parameterization, the or- 
bital phase for a quasi-circular inspiral becomes 



4>{F) 



dE fdE 



dbj \ dt 



jduj 



4>gk{F) 



:A{2TTmF) 



a/3 
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where F and lo ~ inF are the linear and angular or- 
bital frequency, (/>gr = -l/(3277)(27rmF)-5/3 ^he GR 
orbital phase and E{io) ~ — (/i/2)(ma;)^/'^ is the binary's 
binding energy to Newtonian order. Recall here that 
m = mi -|- 7712 is the total mass of the binary, while 
fx = 77117712/777 is the reduced mass and "q = ^/m is the 
symmetric mass ratio. Equation (jl44p is not valid when 
a = 5 (a 2.5PN correction), as then the integrand be- 
comes proportional to uj~^ , which leads to a log term. 



Before we compute the Fourier phase, we must first de- 
fine to J the time at which the stationary phase condition 
is satisfied F{to) = //2, where / is the GW frequency. 
This condition can be solved to yield 



0,GR 



1 
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A{TTmfy'^ 
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where to.GR is the GR t^. Again, this expression is not 
valid at a = 8, because once more the correction to to(/) 
would be a log term. 

With this at hand, we can now compute the Fourier 
phase in the SPA: 



*GW=20(io)-27r/to 



GR 
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where *gr = (3/128)(7r7W/)-^/3, and where M = 
rf/^m is the chirp mass. Again, these expressions are 
not valid when a = 5 or a = 8, for the reasons described 
above. 

The corrections to the GW phase found here map 
directly to the parameterized post-Einsteinian (ppE) 
framework [3l|. In that framework, one postulates that 
modified gravity theories affect the Fourier phase of the 
GW response function in the SPA via 



^ GW 



*GR + /3ppE {T^MfT 



(147) 



where (/3ppE,6ppE) are ppE parameters. We see that this 
is identical to the corrections introduced by a change in 
the energy flux, with the mapping 



/3p 
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16 (a-5)(a-8) 
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-a/ 5 



This is not surprising, as the ppE framework was in part 
motivated by studying power-law (in velocity) modifica- 
tions to the energy flux and the binding energy [3l[. 

We have then found that a large number of energy 
flux corrections associated with extra gravitational and 
scalar fleld emissions can be mapped to the ppE frame- 
work. In the even parity case, the leading-order fre- 
quency exponent 6ppB = —7/3, while in the odd-parity 
case 6ppE = —1/3, unless the binary is non-spinning in 
which case 6ppE = +3. 

The results found in this paper could help in the gen- 
eralization of the ppE framework to more generic quasi- 
circular inspirals. The original framework considered 
only non-spinning, equal mass inspirals, while recently 
Cornish et al. [3^ generalized it to non-spinning, unequal 
mass systems through A — >■ Ar]'^. In this paper we have 
found that A does not only depend on a simple power law 
of 77, but also on the mass difference 5m /m = \/\ — 477 
and on combinations of the spins. For single detections, 
however, such a generalization is not needed as one only 
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Sector A a 

Even-Parity leC^j,^ -2 

Odd-P, Spins (^), ^(4^ [A^ + 2 ((A ■ vi2 f)J +4 

Odd-P, Spins (ft), iC4i(SI5^ (2fi^/ - 3n%^)}^ +4 

Odd-P, No Spin iC4^ +14 



TABLE II: Coefficients of the relative energy flux. 



measures a single number, /?ppE, and one cannot extract 
the dependencies on 77, Sm/m, and the spins. 

Although we currently lack any GW detections, we 
can still estimate the projected constraints that such de- 
tections would place on quadratic gravity. According to 
Table |lll the even-parity sector leads to the strongest 
deviations from GR, since a is the most negative. 
Therefore, we consider EDGE theory, (ai, a2, as, /?) = 
(1, —4, 1, aEDQB)aEDGB, as a simple sub-case of the even- 
parity sector. Let us first imagine that we have detected 
a GW with Ad. LIGO and signal-to-noisc ratio (SNR) 
of 20 that is consistent with GR and that originates 
from a non-spinning BH binary with masses (mi, 1112) = 
(6, 12)M0. Given such a detection, Cornish et al. |32| 
estimated the projected bound |/3ppE| ^ 5 x lO"** for 
^ppE = — |, which implies |q;edgb|^^^ ^5 4 x 10''' cm. 
Let us now assume that we have detected a GW with 
LISA classic with and SNR of 879 and still consistent 
with GR, but that originates from a non-spinning BH bi- 
nary with masses (mi, m2) = (10®, 3 x lOflM© at z = 1. 
Given such a detection, Cornish et al. [32] estimated a 
bound on |/3ppE| < 10"® for the same value of 6ppE as 

before, which leads to ag^Qg < 10^° cm. In both cases, 
notice that these projected bounds arc consistent with 
the small-coupling requirement ^ 1; i.e. saturating 
the projected Ad. LIGO and LISA constraints we have 
CacI. LIGO ^ 3 X 10"^ and Clisa ~ 10"'"' for those par- 
ticular binary systems, which is clearly much less than 
unity. 

Comparing these results with the current constraint 
obtained by the Cassini satellite, |q;edgb|^^^ < 1.3 x 
IQi^cm HI, we see that Ad. LIGO and LISA could 
constrain q;edgb much more strongly. Unfortunately, it 
seems difficult to put constraints on EDGB with binary 
pulsar observations, since NSs have no scalar monopole 
charge in this theory. We emphasize again that this is 
opposite to the expectation from scalar-tensor theories, 
in which NSs have scalar monopole charges while BHs do 
not. Finally, one cannot estimate the bounds one could 
place on dynamical CS gravity, since one would have to 
properly account for modifications to the conservative 
equations of motion, which we have not calculated here. 



VIII. CONCLUSIONS AND DISCUSSIONS 

Wc have studied the binary inspiral problem in a wide 
class of quadratic gravity theories in the slow-motion, 
weak-gravity regime. The structure of a compact object 
in such theories affects the exterior scalar field sourced 
by the object. Despite this, we can model a compact 
object by an effective scalar field source characterized 
by its scalar monopole and dipole moments. The scalar 
monopole charge is enhanced inversely proportional to 
the mass of the object, while the dipole charge is inde- 
pendent of the mass for a fixed dimensionless spin pa- 
rameter. With this effective source, we then derived and 
solved the modified field equations for the scalar field and 
metric deformation. 

We find that the scalar field generically emits dipole 
radiation in the even-parity sector, and quadrupole radi- 
ation in the odd-parity sector. Such radiation affects the 
rate of change of the binary energy at relative — IPN or- 
der in the even-parity case and relative 2PN order in the 
odd-parity case. The quadrupole contribution depends 
quadratically on the BH spins, and thus it is suppressed 
for non-spinning binaries. In that case, the odd-parity 
contribution becomes of relative 7PN order, as found nu- 
merically in [65| . We have found excellent agreement 
between their numerical results and our analytical calcu- 
lations. 

We have also calculated the metric perturbation in the 
FZ and its associated energy fiux. In the even-parity 
sector, the dominant metric contribution leads to a OPN 
relative correction in the energy flux, which is smaller 
than the -IPN correction induced by scalar dipolar radi- 
ation. In the odd-parity sector and for spinning BHs, the 
metric perturbation leads to a 2PN modification to the 
energy fiux, which is of the same order as that induced 
by quadrupolar scalar radiation. In the odd-parity sec- 
tor and for non-spinning BHs, we expect the energy fiux 
correction due to the metric deformation is suppressed to 
at least of 6PN order, as found by Pani et al. (65| . 

Whether these corrections can be measured or con- 
strained depends on whether they are degenerate with 
GR terms in the physical observable, i.e. the waveform. 
A — IPN effect cannot be degenerate, as there are no such 
terms predicted in GR. A 2PN effect, however, could be 
degenerate with a spin-spin interaction for quasi-circular 
inspirals with aligned or counter-aligned spin compo- 
nents. That is, a renormalization of the spin magnitudes 
of both bodies can eliminate this 2PN effect, assuming 
one truncates the waveform at that order. If higher-order 
PN waveforms are used, or if the orbit is more generic 
(i.e. if there is precession or eccentricity), then this de- 
generacy can be broken. 

We also calculated the effects of such energy flux mod- 
ifications on the gravitational waveform. The waveform 
phase depends sensitively on the rate of change of the 
orbital frequency, which in turn is governed by the rate 
of change of energy. We calculated the corrections that 
would be induced in the waveform and mapped them 



22 



to the pp E framework. We then used a recent ppE 
study [32| to estimate the constraints that Ad. LIGO and 
LISA could potentiahy place on quadratic gravity theo- 
ries. Given a GW detection, we found that the magnitude 
of the new length scale introduced by quadratic gravity 
theories (associated with a ratio of their coupling con- 
stants) could constrain at a level controlled by the small- 
est length-scale probed in the inspirals, i.e. the size of the 
smallest compact object's event horizon or surface. The 
best projected bounds achievable with Ad. LIGO will 
thus come from stellar-mass BH or NS inspirals, while 
LISA will benefit the most from EMRIs. Since NSs have 
no scalar monopole charge in EDGB theory, this theory 
cannot be constrained from binary pulsar observations. 
This property is diametrically opposite to scalar-tensor 
theories where BHs have no hair. 

There are several possible avenues for future work. 
Since we here mainly considered corrections due to the 
dissipative sector of the theory, one possibility is to cal- 
culate the non-dissipative corrections that would modify 
the binding energy (here, we mean both gravitational 
and scalar binding energy) and the equations of motion. 
There are two effects that should be accounted for: new 
scalar-scalar forces and metric deformations. Let us con- 
sider the former first. In the even parity case, compact 
objects have an associated scalar monopole charge, and 
thus, there is an additional scalar force with a 1/r poten- 
tial that should lead to a relative OPN non-dissipative 
correction. Similarly, in the odd-parity case, a spin- 
ning compact body possesses a current dipole charge, and 
hence, dipole-dipole interactions should arise. Since the 
dipole potential is proportional to while the dipole 

charge couples to the first derivative of the potential, the 
binding energy and the equations of motion should be 
corrected at relative 2PN order. 

Another non-dissipative modification is induced by de- 
formations of the background metric tensor. In the even- 
parity sector, such corrections enter at relative OPN or- 
der, as found by Yunes and Stein [28j . In the odd-parity 
sector, there is no metric deformation for isolated non- 
spinning BHs, but for spinning ones there is a correction 
proportional to to the {t,i) components Q, which 
then leads to a 4.5PN correction in the equations of mo- 
tion when we consider boosted BHs. This then implies 
the following: (i) in the even-parity case, the conservative 
corrections to the equations of motion do not affect the 
leading-order modification to the waveforms, since this 
is dominated by the — IPN scalar radiation effect; (ii) 
in the odd-parity case, the conservative corrections from 
the metric deformation can be neglected, but those due 
to the scalar-scalar force will contribute at the same order 
as the effect calculated here. A complete analysis of the 
waveform observable would thus require the calculation 
of such a scalar-scalar, conservative effect. 

Another possibility could be to study modified 
quadratic gravity in the context of BH perturbation the- 
ory. This would be a tremendous effort that would have 
to be split into separate parts. First, one would have 



to find an analytic, strong-field solution for arbitrarily- 
fast rotating BHs in quadratic gravity. This has only 
been found in the slow-rotation limit both in the even- 
parity [m and odd-parity sectors @. Once this is ac- 
complished, one would have to study the evolution of 
metric perturbations away from this solution. Such evo- 
lution equations would have to be decoupled in terms 
of some master function to derive Teukolsky-like master 
equations. Finally, with these equations at hand, one 
would have to solve them numerically, when the pertur- 
bations are sourced by a small object in a tight orbit. 
Such an analysis would be interesting because one would 
be able to derive not only the corrections to the energy 
fiux carried out to infinity, but also that which is ab- 
sorbed by the BH horizons and which we ignored in this 
paper. 

A final follow-up would be to study how NS solutions 
are modified in quadratic gravity (sol . and how the en- 
ergy flux from NS binaries is modified. This could then 
lead to direct constraints on quadratic gravity theories 
from double binary pulsar observations. Such constraints 
could be stronger , relative to current Solar System con- 
straints, as they could potentially provide constraints of 
roughly the order of magnitude of the NS radius. Of 
course, in the case of EDGB theory or dynamical CS 
gravity, these constraints might not be stronger as NSs 
have no scalar monopole charge in such theories. 
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Appendix A: The Balding of Neutron Stars in 
EDGE Gravity 

In this appendix, we consider the scalar field equation 
in EDGB gravity for isolated NSs. Integrating the evo- 
lution equation, we find 



(Al) 



where we have defined the Gauss-Bonnet invariant 
7^^B = - 'iR^,yR^"' + Rf.^saR^""^" ■ Since the Gauss- 
Bonnet combination is a topological invariant, the right- 
hand side identically vanishes for any simply-connected, 
asymptotically flat spacetime. Moreover, since we are 
considering isolated NSs, these must be stationary, and 
so the time integration can be removed. 

With all of this and using Stokes' theorem, Eq. (jAl|) 
becomes 



^{d,d)n'dS 



-g{drd)dS = , (A2) 



where is the radial unit vector and the integral is per- 
formed over the 2-sphere at spatial infinity. Notice that 
^/—g ^ r^, while the scalar field must decay at infinity 
for it to have a finite energy. 

Equation (jA2p does not vanish at spatial infinity for 
all scalar field solutions, i.e. if we model i? = z?„/r 
with i9„ a constant, then Eq. (jA2p leads to the unique 
solution -dn = 0. This is a physicists's proof that the 
EDGB scalar field cannot have scalar monopole charge 
for a spherically symmetric NS. Similarly, one can show 
that NSs cannot have scalar monopole charge in dynam- 
ical CS gravity; the proof laid out above carries through 
with the replacement TZq^ — >■ *RR, since *RR is also a 
topological invariant. 



Appendix B: Integration techniques 

In this appendix, we provide some useful integration 
techniques. When computing near-zone integrals, we are 
faced many times with integrals of the form 



d'x^-^ 
rir2 



(Bl) 



When the point-particle approximation is valid, such 
near-zone integrals can be Hadamard regularized by 
keeping only the finite part. Let us then define [68l | 
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to be evaluated in the near-zone and where FP b=o stands 
for the finite part operator (in the limit B 0) and \x\ 



is an analytic continuation factor [68|. The solution to 
this integral is 



1 + 



tE4 

q=0 



The first few Y^^^ are simply 



Yo=Y = b, Y, = -{x\+xl) 

^{ijk) — ^ y^i + Xi X2 + Xi X2 + X2 J 
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(B4) 
(B5) 

(B6) 



The solution to the Y^^,) integral can also be derived 
by using certain Poisson integral identities (63| : 



where we have defined 

Pif) ^ ^ / 
and the boundary term is 
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(B8) 

n^dn' . 



(B9) 

As usual, we retain only those terms that are independent 
of the boundary TZ. 

Finally, there is yet another type of integral that com- 
monly appears in near-zone integration: 



r dV 

\x' - xi\\x' - ; 
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I \iO — xi\\ui; — X2\\x' — X 

Let us then define the so-called triangle potential (69j 

g{Xi,X2,X3) = — — — — -. 

47r \x' - xi\\x' - X2\\x' - X3\ 

(Bll) 

It is a bit of a miracle that the above integral has the 
closed-form solution Q{xa,xb,Xc) = 1 — lnA{ABC), 
with A{ABC) = \xa - xb\ + \xb - xc\ + \xc - xa\- 

One can show that the triangle potential satisfies a set 
of relations, including (69| 
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and more generally 



be rewritten as 



(B13) 



where Q{ABC) = Q{xa-,xb,xc)- 



Appendix C: Odd-Parity, Non-Spinning, 
Regularized Contribution in the Metric Correction 



We consider here the odd-parity sector for non- 
spinning binaries, where, for the scalar field, the 
magnetic-type dipole moment vanishes, = 0, since 
Xa = 0. For the regularized contribution, we only need 
to consider the cross-interaction terms since the isolated 
non-spinning BH solution in the odd-parity case is sim- 
ply the Schwarzschild metric. The k']^^ source term gives 
the largest contribution and one is then left only with the 
pseudo-scalar generated by interaction terms, as given in 
Eq. (HI]). 



The metric deformation is given by Eq. (|95p , the m = 
piece of which can be split as in Eqs. (|110p - (|113p . Before 
tackling each of these terms separately, let us point out 
that many of them identically vanish. For example, one 
of the contribution in Eq. (|110p is proportional to 
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27Tmiejki lim d^gl^^dlpY(xi,X2) 
27Tm2e,kid^^'J,dlf^Y{x,,X2) - . (CI) 



It is critical in this calculation and in the calculations 
that follow to replace the derivatives by particles 
derivatives, i.e. derivatives with respect to x\ and 

Let us then tackle the first contribution to the dissipa- 
tive metric deformation. Equations pi0p - (|113p can then 



,(1) 



: 20487r 



P r 
+ j)l -f (1 ^ 2) , 



u(2) Anna "4 "l?m2 



,(3) 



+ hi\]n + h;\jn + j) +(1^2) 

Anna "4 "^?m2 



(C2) 



(C3) 



(4) 



+(lo2), (C4) 

20487r— vin 

p r 

hijn + hijn + o j) + (1 o 2) , 



(C5) 



where we have defined 
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and 



ji%,c = iimai^)ai^)a(?)e(ABC). 



(C7) 



with A,B,C denoting the multi-index lists. We provide 
a more detailed discussion of J tensors in Appendix |D] 
One can then show through explicit computation that the 



two terms combine to give Inj 4 

0, hijn + I&ijn = 0, and I-!ijyi 
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0. Therefore 



at leading order. 



Let us now look at contributions that are smaller by 
0{v). Such a correction can arise from two different 
terms: (i) the 0{v) correction to the source term with 
TO = in the sum of Eq. ([Ml) , or (ii) the 0{v^) correction 
to the source term with m = 1 in the sum of Eq. (|95l) . 
For case (i), the next-order terms consist of two time 
derivatives and one factor of h^i (or three time deriva- 
tives and one factor of hij)^ which when combined are 
OiiP') smaller than the 0{v'^) contribution shown to van- 
ish previously. Also, the next-order terms in the FN met- 
ric appears at 0{v'^) higher relative to the leading-order 
terms. Finally, i?"^ in Eq. ((55|) expanded as in Eq. (2.27) 
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of [63 with m = 1 in the sum, gives an 0{v) relative 
contribution to d^'i), but exphcit calculation shows that 



d_ 

dt 

d_ 

dt 

= 16-mim2eyfc^ 



qNz 8 a 

TT P dt 

8 a 

TT /3 dt 
,0: 



Vl2k I I — ) ( — ) (f'x 



„ ;)(i)a(2) 



Vl2kdl'^dl^^b 



d-^x 

M fi r2 

: 0. 



(C8) 



For case (ii), the resulting ijij contains one vector. The 
correction to the energy flux consists of multiplied by 
hjj^ and averaged over a 2-spherc. However, since the 
leading contribution in /ijy^ contains even numbers of 
vectors, the correction only contains angular integrals of 
odd numbers of 's which vanish exactly upon integra- 
tion. 

Since there is no 0{G^,v) relative contribution to 
dti)ij, the first, non- vanishing contribution must be at 
least 0{v'^) smaller than what we computed in Eqs. (jC2[) - 
(|C5[) . which amoimts to a 7PN correction to the energy 
flux carried by the metric deformation, in the odd-parity, 
non-spinning case. 



Appendix D: Evaluating J tensors 

Recall that the definition of the J tensors is 

Ji^kc - lini d<^^d<^^^dj^^g{ABC) . (Dl) 

' 3— 



Another type of problematic limit is 



lim 



or lim ' 

3^p 



''p3 : 



(D3) 



which does not formally exist, since it depends on the 
path taken as we describe below. Parameterize the path 
that particle 3 takes to the location of particle p by the 
continuously differentiablc path 7(A), with A a parameter 
of path length and A = the location of particle p. There 
are an infinite number of paths one could choose, and 
each can be parameterized in two senses. Taking the 
limit along this path "from below" (i.e. from smaller 
values of A to larger values) yields 



'pS 



-«;(o), 



(D4) 



lim 

3— ^p, 7~ 

where is the tangent vector to the curve 7. Taking the 
limit from above, we find 



lim 

3— ^p, 7"^ 



^3 



<{0) ■ 



(D5) 



The limit depends on the path's tangent at the point of 
particle p, and the direction in which the limit is taken. 
Clearly, the final answer must be unique, which implies 
the limit must vanish. 

A unique prescription to this problem is formalized as 
Hadamard rcgularization [to} . This can be summarized 
as follows. All possible paths are considered, with tan- 
gent vectors v-y. The average is then taken by integrating, 
e.g. 



The limit 3 — >■ p which appears must be taken with care. 
There may be terms proportional to 



lim — 

3^-p rp3 



which have no finite part. In the evaluation of the J 
tensors, only the finite part of the limit is kept. That 
is, a function can be expanded as a Laurent series about 
these points, and the finite part scales as (rp3)° in the 
limit as 3 — > p. 



lim • • • n;^3 ■ 

3— >p 



47r 



.ylJ... 



(D2) The first few such limits, for example, are 



limn;3=0, 



lim n 

3-)-p 



p3 



(D6) 

(D7) 
(D8) 
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